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CHEMICAL  SYNTHESIS  UNDER  HIGH  PRESSURE  SHOCK 
COMPRESSION  LOADING  :  MATERIALS,  MEASUREMENT  AND 

MODELING. 


R.A.Graham 

(Sandia  National  Laboratories,  Albuquerque  ,The  Tome  Group  ,  Tome,  New 

Mexico,  USA) 

It  is  now  widely  recognized  that  synthesis  of  compounds,  alloys  and 
amorphous  materials  can  be  readily  accomplished  with  the  application  of  high 
pressure  stress  pulses  produced  by  intimate  contact  with  detonation  waves  from 
high  explosives  or  high  velocity  impact.  Nevertheless,,  after  30  years  of  research,  it 
is  apparent  that  there  is  no  consensus  as  to  the  appropriate  fundamental  descriptions 
of  the  shock-induced  chemical  synthesis  processes.  Indeed,  fundamentally  different 
and  contradictors'  paradigms  are  invoked  by  various  investigators.  The  basic  choices 
for  scientific  definition  of  the  processes  are  paradigms  based  on  either  mechanical- 
and  thermodynamic-equilibrium  or  nonequilibrium  mechanochemical  processes.  The 
present  paper  summarizes  the  status  of  shock-induced  chemical  synthesis  science 
and  indicates  directions  required  to  develop  realistic  scientific  models  of  powder 
mixtures  in  the  highly  porous  (about  50%)  state. 

Based  on  the  newly  available  nanosecond,  time-resolved  stress  and  stress- 
rate  measurement  capabilities  as  well  as  combined  PVDF-VISAR  techniques,  there 
IS  now  specific  quantitative  date  which  to  develop  scientific  models  of  the  shock- 
compression  processes  resulting  in  chemical  synthesis.  The  data  on  titanium-silicon, 
aluminum-hematite,  and  inert  mixtures  provide  rich  detail  of  the  details  of  the 
mechanochemical  processes.  It  will  be  necessary-  for  fully  qualified  models  to 
incorporate  coordinated  studies  of  materials  science,  nanosecond  time-resolved  wave 
profile  measurements  and  mechanochemical  models  which  are  based  on  qualitative 
and  quantitative  data  on  a  range  of  materials  in  various  porosities  and 
morphologies.  It  appears  that  the  field  of  shock-induced  solid  state  chemistry'  is  a 
watershed  from  an  early  exploratory  phase  to  a  more  quantitative  scientific  phase 
from  which  detailed  models  can  be  developed. 
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DYNAMIC-STATIC  COMPRESSION  OF  INORGANIC  MATERIALS 


S.S.  Batsanov 

(Center  for  High  Dynamic  Pressures,  Mendeleevo,  RUSSIA) 

Introduction 

40  yeas  ago  Yu.N.Rabinin  [1]  for  the  first  time  had  preserved  and  studied  the 
shocked  chemical  compound  and  by  this  opened  a  new  chapter  in  science,  namely  the 
shock  induced  chemistry. 

Development  of  this  field  permitted  to  carty  out  constructions  of  the  different  recovery' 
ampoules,  to  realize  the  physical-chemical  transformations,  to  obtain  the  new  phases  and 
materials.  Results  of  these  investigations  are  reviewed  in  books  [2-5], 

In  the  shock  recovery'  experiments  the  temperature  may  be  varied  from  a  few  degrees  to 
several  thousand  degrees  and  the  pressure  may  be  a  few  kbar  or  a  few  Mbar.  However, 
high  temperature  generated  by  the  residual  heat  lasts  much  longer  than  high  pressure  and 
leads  to  the  annealing  of  metastable  phases,  thus  strongly  devaluating  the  advantage  of 
shock  chemistiy'. 

The  many  methods  of  countering  the  negative  influence  of  the  temperature  factor 
reduce  to  cooling  the  experimental  devices  before,  during  or  after  shock  compression.  As 
well  as  trivial  cooling  of  the  recoveiy  ampoules  by  cryogenic  liquid  prior  to  explosion  or 
the  use  of  explosion  in  water  of  liquid  nitrogen,  other  common  approaches  are  adding 
water  to  the  initial  powder,  pressing  it  in  a  copper  matrix,  applying  a  thin  layer  between 
massive  metal  plates,  or  else  the  detonation  of  a  mixed  charge,  where  temperature 
reduction  accompanies  scattering  of  the  detonation  products. 

Another  way  of  reducing  the  annealing  effects  of  a  heat  is  to  extend  the  duration 
of  the  dy'namic  pressure  by  increasing  the  charge  or  the  density  of  the  surrounding  media 
(water,  ground,  metal  shells).  However,  these  techniques  extend  the  duration  of  the 
pressure  by  only  1-2  orders  of  magnitude  and  thus  have  no  significant  influence  on  the 
temperature  parameters  of  the  unloading  body  and  hence  on  its  structure  and  properties. 

In  our  laboratory'  was  cany'  out  a  novel  approach  to  this  task  -  the  method  of 
dynamic-static  compression  (DSC). 

Principle  of  the  DSC-method 

If  a  strong  recovery'  ampoule  contains,  besides  a  sample  to  be  compressed  a 
condensed  substance  (for  example,  water)  able  to  evaporate  under  action  of  the  residual 
temperature,  then  after  an  impact  loading  such  an  ampoule  becomes  an  autoclave,  the 
dynamic  pressure  being  replaced  by  a  static  one  [6]. 

However,  opening  of  siich  an  ampoule  (still  containing  high  pressure)  presents  a 
very'  complicated  technical  problem  due  to  its  brisance.  Therefore  we  proposed  [7]  to  use  a 
solid  working  body  undergoing  during  unloading  a  polymorphous  transformation  or 
melting  with  an  increase  of  volume  by  10-20%,  w'hich  is  not  enough  to  provide  a  big 
brisant  effect  (Fig.l).  Recently,  it  has  been  concluded  that  a  high  residual  pressure  may  be 
produced  in  the  absence  of  phase  transformations  in  the  working  body, 
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simply  as  a  result  of  the  difference  in  compressibility  of  the  ampoule  shell  and  its  contents 
[8]  (Fig.2). 

Obviously  a  magnitude  of  the  residual  pressure  is  determined  by  strength  and 
rigidity'  of  the  ampoule  shell.  Partially  one  can  obtain  the  static  pressure  about  15  to  20 
kbar  in  the  strong  cylindrical  ampoule. 

Special  case  is  using  the  explosive  substance,  for  example  RDX,  as  a  working 
body.  If  to  add  some  percents  of  high  explosives  to  a  powder  to  be  shock  compressed,  then 
aher  an  explosion  of  the  outside  charge  we  can  obtain  as  a  result  of  the  inner  detonation  a 
solid  foam.  RDX  can  be  added  as  a  homogenic  mixture  or  in  form  threads,  spirals,  plats. 
In  these  cases  one  can  obtain  solid  foams  with  canals  and  gaps  of  the  special  form. 

As  long  time  one  can  keep  the  residual  pressure?  If  the  working  body  is  a  solid 
the  static  pressure  may  be  conserved  during  few  months.  In  cases  of  liquid  and  gas 
\^’orking  bodies  very'  difficulty  to  keep  the  residual  pressure  under  high  temperature,  but 
after  cooling  the  recovery'  ampoule  to  the  room  temperature  the  static  pressure  reduces 
usualh  to  few  kbars  due  to  a  condensation  some  gases  or  solidification  some  liquid. 

Thermodynamic  parameters  of  DSC-amooules 

DSC>method  was  realized  in  the  cy'lindrical  and  plane  recovery'  ampoules. 
Cylindrical  ampoules  with  a  spherical  end  are  used  from  the  very'  beginning  of 
experiments,  to  eliminate  stress  concentrators.  The  ratio  of  the  diameter  (50  mm)  to  the 
height  (70  mm)  is  0.7,  which  is  such  that  tlie  strength  if  the  cy'linder  approximate  the 
characteristics  of  a  sphere.  The  ampoules  are  exploded  on  a  lathe  with  a  steel  tube  with  an 
aperture  for  cutting,  pre\  enting  the  loss  of  parts  of  the  ampoule  on  thinning  of  the  shell. 
The  ampoule  breaks,  as  a  rule,  when  the  wall  thickness  reaches  1-3  mm;  the  fracture 
thickness  is  a  good  measure  of  the  residual  pressure. 

Increasing  the  static  pressure  in  DSC-ampoule  abo\e  the  le\el  of  15  kbar  is  a 
special  problem,  which  is  may  be  approached  by  several  ways.  A  first  method  is  to 
produce  a  tw'o-chamber  ampoule  in  which  a  static  pressure  P=15  kbar  acts  between  the 
inner  and  outer  shells  and  the  total  pressure  ~  30  kbar  acts  on  the  investigated  material 
(Fig.  3).  The  second  approach  is  to  use  a  working  body  that  contains  several  components 
w  ith  increasing  phase-transition  pressure  or  compressibility'  w'hich  are  either  pressed  into 
a  homogeneous  mixture  or  else  are  placed  individually  in  concentric  layers  within  the 
ampoule  so  that*are  successively  included  in  the  creation  of  the  static  pressure  in  the 
course  of  dynamic  unloading  and,  remaining  inert,  apply  in  the  additional  shell.  Both 
approaches  w'ere  investigated  in  our  laboratory,  however  w'e  think  that  the  maximal 
perspeclix  e  has  the  magnetic  method. 

It  is  know  n  that  the  metal  shell  under  the  magnetic  field  of  10'^  ersted  may  be 
compressed  for  several  tens  microsecond  to  pressure  ca.50  kbar.  The  source  of  magnetic 
field  may  be  the  battery  of  ccndensators  or  the  magneto-cumulatixe  generator.  The 
simultaneous  action  of  shock  wax'es  and  a  magnetic  field  prolongates  the  dynamic 
pressure  till  several  tens  or  even  hundred  microsecond.  We  carried  out  the  simulation  of 
this  process  and  began  to  prepare  the  physical  experiment. 
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Studying  shock-wave  parameters  in  the  recovety  ampoules  is  the  very’  complicate 
task  and  now  this  work  is  in  progress.  The  preliminary'  results  are  presented  on  Fig.4.  As 
one  can  see  the  shock-wave  velocity  approaches  the  detonation  velocity  only  at  the  very' 
bottom  of  the  DSC-ampoule.  The  sample  to  be  compressed  we  put  just  at  this  place. 

To  clarify  the  picture  of  shock-wa\'e  propagation  in  such  ampoules  we  carried  out  the 
experiments  where  the  initial  powder  of  the  working  body  was  intercalated  with  se\’eral 
sheets  of  steel  or  copper  foil  perpendicular  to  the  ampoule  axis.  The  diameter  of  the  hole 
punched  in  such  a  foil  during  shock  compression  is  known  to  correlate  with  the  Mach 
stream  dimensions  [9]. 

In  DSC-experiment  a  hole  in  the  foil  appeared  for  the  first  time  at  a  distance 
from  the  top  of  the  ampoule  equal  to  its  inner  diameter,  Down  from  this  level  the  liole 
diameter  increases  and  then  decreases  (Fig.5),  indicating  the  instability  of  shock  waves' 
interaction  in  the  axial  part  of  the  ampoule. 

Taking  in  consideration  the  formation  of  holes  in  the  axial  part  we  carried  out 
measurement  of  shock  pressures  by  a  manganin  gauge  in  the  form  of  a  ring  (diameter  10 
mm)  which  show  that  P=190  kbar  at  a  distance  of  5  mm  from  the  ampoule  axis  was  used 
an  explosive  with  D=7.6  km/s. 

The  residual  temperature  in  DSC-ampoules  was  measured  calorimetrically.  The 
ampoules  were  plunged  into  a  water  calorimetr  in  45-100  s  after  an  explosion,  then  after 
cooling  the  same  ampoules  were  heated  in  an  electric  furnace  to  400®C  and  placed  in  an 
explosition  chamber  ( in  reproduce  the  cooling  conditions),  where  a  cooling  curves  were 
recommended.  By  this  procedure,  the  residual  temperature  for  zero  delay  after  explosition 
was  determined  by  an  extrapolation  to  an  accuracy  of  5-7  %  .  Calorimetric  experiments 
were  conducted  for  the  actual  ampoules,  for  steel  dummies,  and  for  ampoules  with  w  eakly 
reinforced  bottom  that  were  blown  away  (with  the  filler)  in  the  course  of  explosition.  This 
permits  to  determine  the  shell  temperature  and  to  calculate  the  filler  temperature:  Tres 
=900T  for  D=  7.2  and  1200  °C  for  D=7.8  km/s.  If  a  hole  of  3-mm-diamcter  is  made  in 
the  ampoule  bottom,  the  compressed  material  escapes  from  the  ampoule  (into  the  v  olume 
of  a  coaxial  cylinder  of  the  same  cross  section  in  the  course  of  unloading.  The  residual 
temperature  of  this  ampoule  is  measured  by  comparison  with  temperature  of  the  closed 
ampoule  to  estimate  Tres  of  the  axial  region  (the  Mach  stream),  which  was  found  to  be 
1500-2000  T. 

In  our  laboratory'  were  constructed  and  used  the  plane  recovery'  ampoules  for 
DSC-method.  Thennodynamic  parameters  in  these  ampoules  may  be  calculated  with  a 
good  accuracy. 


Results  of  DSC  of  inorganic  materials 

The  described  method  was  used  for  the  inv  estigation  of  the  BN  phase  transition 
in  the  plane  and  cylindrical  ampoules. 

The  plan  e.xperiment  was  realized  by  using  the  gas  gun  (Fig.  6)  [10].  The 
thickness  of  sample  is  equall  mm,  one  of  projectile  is  equal  5  mm.  the  weight  of 
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projectile  is  equal  150  g.  Time  of  impact  is  ca.  2.5  micro-second.  Conditions  and  results 
of  ballistic  e.xperiments  are  given  in  Table: 


Velocity  of 
projectile, 

Pressure,  kbar 

Yield  of  w- 
BN,  % 

m  /s 

Shock 

Residual 

673 

156 

0 

1  5-20 

678 

158 

0 

15-20 

685 

160 

10 

80-90 

As  one  can  see,  the  residual  pressure  increases  the  yield  of  the  dense  phase  of  BN  few 
times. 

Using  the  cylindrical  ampoule  and  a  melt  of  trotyl/hexogen  (40/60)  permitted  to 
realize  the  complete  phase  transition  h-BN  in  w-BN  and  form  a  strong  monolithic 
sample  of  specified  dimensions  [11].  The  material  produced  has  a  microhardness  of 
5000-7000  Vickers  units  in  the  axial  region  (in  the  form  of  an  especially  hard  spot 
extending  over  up  to  15%  of  the  internal  diameter  of  the  ampoule),  which  falls  to  1000- 
1500  units  at  the  periphery  of  the  cylindrical  sample;  Fig. 7.  The  sample  density  is  3.35 
g/cm^  ,  and  the  pressure  of  crushing  is  110  kg/mm^  .  X-ray  data  show  that  the  sample 
hardness  varies  antibatically  with  variation  in  the  regions  of  coherent  scattering  and  the 
domain  size  is  ca,  100  A  for  the  hardest  part  of  sample. 

DSC  of  diamond  powder  also  leads  to  the  formation  of  a  hard  briquette  with  a 
microhardness  of  up  section.  Approximately  the  same  in  samples  obtained  by  the  DSC  of 
wlirtzitlike  and  cubic  form  BN  powder.  A  microhardness  all  these  materials  are  ca. 
5000-6000  Vickers  units  in  the  axial  region,  densities  ca.  3.3  g/cm^  .  DSC-method  was 
used  with  a  success  for  compacting  of  SiC,  B4C,  B,  and  TiB2  powder. 

DSC  the  mixtures  of  powders  of  Si02,  GeOs,  SnOs  and  10-35%  RDX  lead  to  a 
formation  of  solid  foams  of  these  materials .  X-ray  analysis  showed,  that  the  Si02  -foam 
has  the  quartz  structure,  Ge02  -tetragonal  structure,  and  Sn02  is  cassiterite.  The 
density  of  the  Si02  -  foam  varies  from  0.6  to  1.4  g/cm^  ,  densities  of  the  Ge02  and 
Sn02  -foams  are  2.6  and  2.9  g/cm^  respectively. 

An  estimation  of  temperatures  in  these  experiments  leads  to  magnitude  1200  to 
1500  °C.  This  temperature  is  enough  for  the  contact  melting  and  bonding  of  particles  and 
a  formation  of  solid  foams. 

An  attempt  to  farm  the  foams  from  Mo,  Si  and  MoSi2  by  the  same  conditions  was 
unsuccessfully.  The  reaction  consists  in  the  high  temperatures  of  melting  of  these 
materials.  When  we  used  for  DSC  the  mixture  of  composition  Mo+  10%  RDX,  the  foam 
was  syntheses  very  easy.  Probably  the  additional  temperature  from  the  exothermic 
reaction  of  the  MoSi2 
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formation  was  enough  for  melting  of  a  product  and  bonding  particles,  This  fact  shows  that  the 
chemical  reaction  is  undergoing  during  the  shock  compression. 

We  carried  out  also  DSC-  experiments  with  mixtures  of  the  SiOt  powder  and  small  additions 
of  the  different  metals  (  Cu,  Nu  Pt  )  and  oxides  (EioOj,  Ti02)  pow'ders,  which  may  be 
catalysts.  Taking  into  consideration  that  by  Ii)SC -method  we  can  obtain  porous  and  additions 
ot  different  forms  and  concentrations  in  the  any  solid  foams,  this  method  may  ha\'e  a  practical 
applications. 
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Fig„  Estimation  by  the  foil  method  of 
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Fig*6  Schema  of  ballistic  experiments  with  the  plane  DSC~ampoule 
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COMPUTER  AIDED  STUDY  OF  REACTION-ASSISTED 
POWER  MIXTURE  SHOCK  COMPACTION  AT  MESO-SCALE. 

NEW  COMPUTATIONAL  TECHNIQUE. 

S.G.Psakhie,  E.  V.Shilko  ,  A.  Yu.SmoIin,  A. I. Dmitriev,  S.  Yu. Korostelev 
(Institute  of  Strength  Physics  and  Materials  Sci.,  Rlis.  Acad.of  Sci.. Tomsk,  Russia) 

1.  MCA  method  general  formalism 

The  method  of  Movable  Cellular  Automata  (MCA)  [1-3]  is  developed  on  the  base 
of  rigid  body  dynamics.  At  the  same  time  it  includes  all  the  advantages  of  conventional 
cellular  automata  approach.  This  has  become  possible  due  to  introduction  of  the  new 
definition:  -  *'the  state  of pair  of  automata**  (“relation”  of  interacting  pair  of  automata)  in 
addition  to  the  conventional  one  of  “the  state  of  a  separate  automaton”.  Note  that 
introduction  of  this  definition  allowed  to  go  from  the  net  concept  to  the  concept  of 
neighbors,  that  is  of  principal  importance.  As  a  result  of  this  the  automata  have  ability  to 
change  their  neighbors  by  switching  the  states  (“relations”)  of  the  pairs  of  automata. 

Thus,  within  the  new  method  of  Movable  Cellular  Automata  the  automata  are 
allowed  to  move,  owing  all  the  advantages  of  classical  cellular  automata  approach. 

Within  the  framework  of  the  Movable  Cellular  Automata  method  the  simulated 
material  is  considered  as  an  ensemble  of  discrete  elements  (cellular  automata),  interacting 
by  certain  rules,  relations  and  laws.  Such  automata,  in  turn,  combine  the  elements  of 
heterogeneous  media,  in  particular:  separate  grains  of  polycrystal,  separate  particles  of 
powder  mixture,  etc.  Certainly,  the  size  of  the  automata  are  defined  by  conditions  of 
specific  problem.  It  is  very  important  to  select  the  size  of  automata  correctly  for  adequate 
description  of  the  material. 

Mechanical  state  of  the  simulated  system  is  characterized  by  the  following 
parameters:  position-vectors  of  the  elements  {/?'};  translation  velocities  of  the  elements 

{F' } ;  rotation  angles  of  the  elements  {6']  ,  rotation  velocities  of  the  elements  {tw ' }  . 
Ever^^  automaton  is  characterized  by  the  size  parameter  d ',  mass  w'  and  inertia  moment 

tensor  J  ^ .  Mechanical  behavior  of  the  ensemble  of  the  automata  is  defined  by  “relations” 
of  interacting  pairs  and  the  forces  of  inter-element  interaction.  The  type  of  relations 
defines  tlie  type  of  forces  interacting  between  automata.  In  the  simplest  case  the  relations 
among  the  pair  of  interacting  automata  may  be  divided  in  the  two  types:  linked  (the  pair 
of  the  automata,  among  which  there  are  any  chemical  bonds)  and  unlinked  (there  is  no 
any  chemical  bonds  among  the  automata  of  the  pair).  Note  that  the  unlinked  pairs  of 
automata  may  be  divided  in  its  turn  in  the  two  types:  contacted  (the  pairs  of  interacting 
automata)  and  uncontacted  (there  is  no  any  interaction  among  the  automata  of  the  pair), 
when:  Ir’  =  d’  -  rf  >  0 

Formation  of  initial  structure  is  performed  by  setting  up  certain  relations  to  each 
pair  of  neighboring  elements.  The  change  of  relations  among  the  cellular  automata  is 
defined  by  the  value  of  their  overlapping  h'^ .  calculated  as  h'-’  =  -  rf  ,  where  is  the 
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distance  between  the  centers  of  the  neighboring  elements,  and  rj-'  is  defined  as  rj-’  =  (</ 
'  +  (}■’  )/2,  where  d  ’  is  the  size  parameter  of  the  automaton.  The  pair  of  linked  automata 
remains  to  be  linked  while:  /?''  <  h'-’^ax  The  pair  of  unlinked  elements  remains  to  be 
unlinked  while;  ^  Zf  ■  If  the  two  elements  were  initially  unlinked^  they  become  to  be 
linked,  when;  h-'  <  .  And  vice  versa  -  if  the  two  elements  were  initially  linked.  tlie>^ 

become  to  be  unlinked,  when  their  overlapping  parameter  h  It  must  be  noted 

that,  if  mutual  rotation  of  the  elements  are  taken  into  account,  then  it  is  necessaiy  to 
introduce  the  criterion  of  transition  of  the  pair  from  the  linked  to  unlinked  relations  due 
to  their  relative  rotation.  In  this  case,  if  A  f^\p>  A  r\a.-c  ,  then  the  pair  been  linked 
becomes  unlinked.  Here  A  is  the  distance  between  the  contact  points  of  the  two 
interacting  automata,  calculated  from  the  moment,  when  they  become  linked. 

The  pair  approach  is  used  for  description  of  the  inter-automata  mechanical 

interaction.  is  the  force  acting  from  j-th  automaton  to  Z-th  one.  In  general  case  F'^ 
can  be  represented  as;  F’^  =  ,  where  p’^  is  central  type  of  interaction,  and 

is  a  tangential  type. 

The  central  force  of  pair  interaction  may  be  divided  into  two  components;  elasto- 
plastic  component  (depending  on  degree  of  deformation)  and  viscous  component 
(damping).  The  elasto-plastic  component  reflects  the  property  of  the  material  to  resist  the 

deformation  and  is  defined  by  degree  of  deformation  of  tlie  automaton  F  -  - — 

d'jl  ’ 

where  q'  is  the  distance  from  center  of  automaton  to  contact  point  of  this  automaton  with 
its  neighbor.  In  general  case  the  following  regions  may  be  resolved  when  calculating  the 
force  resisting  to  deformation  (Fig.  1); 


Fig.  1.  Force  resisting  to  deformation. 
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1.  The  region  of  eJaslic  interaction  +  A{^y,l  where  S,  is  the  limit  of 

elasticity  of  the  material  under  compression  (c  -  compression)^  and  As'y,  is  that  one 
under  tension  (/  -  tension). 

2.  Tlie  region  of  plastic  flow  under  tension  \.f!o  +  Ac'y, :  wlicre  f'o^is  the  strength  limit 

of  (he  material  under  tension. 

3.  The  region  of  plastic  flow  under  compression  \^o  -  where  is  the  ultimate 

compression  deformation.  This  parameter  is  defined  by  the  strength  of  the  material 
under  compression  and  its  compressibility.  The  cui^^c  corresponding  to  central  force  is 
labeled  1. 

In  the  region  of  elastic  interaction  the  force  of  resistance  in  the  simplest  case  may 
be  calculated  following  the  linear  dependence: 


(2) 

where  5^^  is  the  contact  area  of  the  automata  /-th  and y-th  ones,  is  the  modulus,  defining 
the  materials  response  to  shape-changing  under  deformation. 

In  the  simplest  case  the  value  of  E  may  be  defined  by  such  a  way.  Consider 
uniaxial  tension  (compression)  of  a  parallelepiped  shape  specimen  along  X-axis.  Assume 
that  it  is  deformed  uniformly.  Then  following  X-component  of  stress-tensor  will  be: 


AA" 

~Y 


1  AK' 


^3  V 


+  K- 


AV 


(3) 


where  AV is  the  volume  change  under  deformation,  G  is  the  shear  modulus,  K  is  the  bulk 
modulus.  Since  deformations  are  uniform  AY=AZ=-c7AX,  and  then  -2cr)/lXX., 

where  cr-  is  Poisson  ratio.  So  the  stress  corresponding  to  resistance  to  shape  change  will 
be 


r.  =  2G  1 


1  ~  2c7^  AX 


And  we  can  define  our  modulus  as  j?'  ^  2G‘ 


1  - 


3 

1  -  2o-' 


X 


In  the  regions  of  plastic  flow  in  the  simplest  linear  case  tlie  resistant  force  may  be 
calculated  as: 


E{Pi){B)  -  -  [E'Ae! Yt  +  E’lZB  -  Bq  - 

where  E’u  is  the  coefficient,  taking  effectively  into  account  the  hardening  of  the  material 
of  /-th  automaton  under  plastic  flow. 

The  viscous  component  describes  the  dissipate  properties  of  the  material  and  is 
defined  by  the  velocity  of  relative  movement  of  the  elements  ^  „  , where 

is  the  coefficient  of  the  force  of  normal  damping.  is  the  normal  component  of 
the  relative  velocity  ofy-th  automaton  in  the  coordinate  s)^stem  centered  at  /-th  automaton. 
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The  langenlia!  force  of  pair  interaction  may  be  divided  into  three  components; 
elasto-plastic  component  (depending  on  degree  of  shear  deformation),  viscous  component 
and  dn,’  friction  force  (only  for  unfinketl  and  contacted  pairs  of  automata). 

Within  this  model  the  force  of  resistance  to  shear  deformation  is  applied  to  the 
linked  automata  only.  The  angle  of  shear  deformation  y  '  of  /-th  automaton  under  its 
interaction  with y-th  automaton  is  defined  by  relative  translation  movement  of  /-th  and y-th 
automata  and  rotations  of  y-th  and  /-th  automata.  Increment  of  shear  deformation  angle  y' 
is  defined  as  follows.  The  total  tangential  displacement  of  /-th  and  y-th  automata  relative 
to  their  contact  point  in  one  time  step  At  will  be 

A/ -  S'q'  -  co^q^)  x  , 

where  x  (F'  -  V'  )jr'^  .  This  displacement  is  considered  as  sum  of  increments 

of  shear  deformation  angles 

A/'^  =  q '  tan( A/ ' )  +  q^  tan(A>'  ) . 

therewith  the  following  equation  is  required;  f'^  {y')-  -f^'  (/  *' )  To  oppose  rotation  of 
the  automata  as  gears  {co'^  =0,  eb'q'  =-d}^q^)  the  couple  forces  is  applied  to  each 

automaton.  Tliis  forces  are  proportional  to  O'  —  0^ . 

The  expression  for  calculation  the  viscous  friction  force  may  be  written  as; 

where  rj^  is  the  viscous  friction  coefficient,  is  the  tangential 

component  of  the  relative  velocity  of y-th  automaton  in  the  coordinate  system  centered  at  /- 
th  automaton. 

Dry  friction  force  flip'-’ )  if  =  0  is  varied  from  zero  to  its  maximum  value 
.  This  force  is  equal  to  tangential  component  of  resultant  force  of  reaction 
defining  by  all  the  forces  acting  on  the  automata.  Here  y?-'  is  vertical  component  of  the 
force  F^.  If  the  value  of  the  tangential  component  of  the  force  F^  exceeds  the  value  of 
p'’p'\  then  and  the  dry  friction  force  reaches  its  maximum  value. 

Resistance  of  the  automata  to  the  change  of  its  volume  is  taken  into  account  by 
analogy  to  the  “embedded  atom”  method.  Volume  dependent  force  Fp’ ,  acting  of  /-th 
automaton  due  to  environment  pressure,  may  be  calculated  as; 

j 

where  the  index y  (/’  =1...V)  numbers  the  automata  being  linked  or  contact  with  the  /-th 
one;  is  the  pressure  of y  -th  neighboring  automaton;  F’  is  the  contact  area  of  /  -th  and 

y  -th  automata;  is  the  unit  vector  normal  to  the  contact  area  and  directed  to  the  y-th. 
The  pressure  P  ^of  y  -th  automaton  is  defined  by  change  of  the  volume  of  automaton,  and 
the  volume  change  may  be  calculated  from  changes  of  the  distances  from  center  of  y  -th 
automaton  to  contact  points  of  this  automaton  with  its  neighbors.  In  Fig.  1  the  cun^e 
corresponding  to  volume-dependence  force  is  labeled  2. 
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Since  in  MCA  method  the  size  of  an  automaton  is  introduced  directly,  then  for 
complete  description  of  the  system  evolution  it  is  necessary  to  add  the  equations  for  the 
moments  to  the  motion  equations  for  the  mass  centers.  Then  the  motion  equations  will  be: 


J 


dt' 
dt^ 


where  =  [p'^  +  x  q'^  is  the  distance  to  contact  point  of  /-th 

automaton  with y-th  one.  The  unit  vector  is  defined  as  =  {R  '  -  R‘  )lr‘-’  ,  where 
is  the  distance  between  the  centers  of  the  automata. 

Using  boundary  conditions  of  different  types  (hard,  elastic,  visco-elastic,  etc.)  it  is 
possible  to  imitate  different  properties  of  surrounding  medium,  containing  the  simulated 
system.  It  is  possible  to  model  different  modes  of  mechanical  loading  (tension, 
compression,  shear  strain,  etc.)  by  setting  up  additional  conditions  at  the  boundaries. 

2.  Chemical  reaction  modeling 

Model  of  virtual  crystal.  In  the  simplest  case  the  cellular  automaton  may  be  of  the 
material  of  the  same  type  (an  initial  compound  of  the  mixture  or  the  product  of  reaction). 
Lets  call  the  state  of  this  element  the  "simple  state"  for  it  is  totally  defined  by 
characteristics  of  the  only  one  compound.  In  general  case  a  cellular  automaton  can  be 
considered  as  the  mixture  of  several  different  substances  and  it's  complicated  state  will  be 
defined  by  the  properties  of  the  all  compounds  and  their  interactions.  To  describe  the 
properties  of  such  a  combined  cellular  automaton  we  follow  the  model  of  virtual  crystal. 
In  this  model  all  the  specific  characteristics  of  an  element  as  well  as  parameters  of  it's 
interaction  with  the  neighbors  are  defined  by  averaging  over  the  all  atoms  of  the 
compounds  making  up  the  element; 


where  bi  and  n,  are  some  certain  unit  characteristics  and  quantity  of  i-th  compound 
correspondingly,  n  is  the  total  number  of  atoms,  N  is  the  number  of  all  possible 
compounds. 

In  this  approach  some  certain  components  of  the  material  composing  an  element 
are  able  to  realize  independently  phase  transformations,  de-solution  into  initial 
components  (in  the  case  of  chemical  compounds)  and  chemical  reactions  among 
themselves.  In  parallel  account  of  these  possibilities,  phase  transformations  will  have  the 
highest  priority,  and  chemical  reactions  will  be  considered  as  the  last  possibility. 

Thus,  movable  cellular  automaton  may  change  its  state  both  because  of  inner 
transformations  and  under  chemical  reactions  with  its  neighbors. 
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Mechanism  of  shock-induced  solid  phase  reactions.  In  general  case  tlie 
mechanism  of  a  solid  phase  reaction  in  a  sliock  wave  is  presented  as  follows.  Only  the 
pairs  of  automata  (of  reagents  A  and  B)  in  the  slates  ’'linked'’  or  “contact”  are  taken  into 
account.  High  compression  of  the  interacting  elements  lakes  place  under  sliock  wave 
propagation,  that  corresponds  to  input  of  elastic  energ)'  into  material.  The  processes  of 
abnormal  high  rate  mass  transfer  and  mixing  the  reagents  are  possible  in  heterogeneous 
media  on  different  scale  levels  (from  atomic  to  micron)  under  similar  conditions. 

The  mechanisms  of  mass  transfer  are  different  and  connected  not  only  with 
intensive  plastic  deformation,  but  also  with  the  grinding  of  the  material  and  further 
mixing  of  the  fragments.  Beside  this,  under  these  conditions  an  intensive  mass  transfer 
may  lake  place  on  the  atomic  level  because  a  portion  of  the  atoms  (especially  in  micro¬ 
damage  zones,  “hot-spots",  etc.)  may  have  the  energ>'  considerably  higher  than  the 
diffusion  activation  threshold.  This  mixture  containing  the  atoms  and  molecules  of  the 
reagents  with  high  kinetic  energies  is  reacting  by  generation  of  the  products. 

In  the  case  of  exothermic  reactions  there  is  characteristic  time  (the  time  of 
transition  to  thermal  equilibrium).  During  this  time  the  released  energy  is  transformed 
into  the  energy  of  mechanical  compression  of  neighboring  regions  and  into  the  thermal 
energy  of  atomic  oscillations.  The  formation  of  the  reaction  product  is  completed  when  in 
the  whole  automata  volume  the  processes  of  mixing,  reacting  in  the  mixture  and 
redistribution  of  the  released  energy  are  finished. 

Reaction  between  the  pair  of  elements.  It  is  possible  to  divide  the  process  of  the 
reaction  of  a  pair  of  the  automata  into  3  main  sub-processes,  ongoing  in  parallel; 

1 .  Mixing  the  reagents  and  formation  stoichiometric  mixture. 

To  estimate  the  mixing  velocity  the  following  equation  may  be  used; 

where  ~  aw  upper  limit  of  penetration  velocity  of  A(B) 

component,  is  some  effective  coefficient,  is  the  elastic  energy  of  compression  of 

A(B)  element  per  1  mole,  is  the  binding  energy  of  the  atom  of  A(B)  component, 

/Jia{b)  is  the  molar  mass. 

2.  Forming  the  product  of  the  mixture. 

To  estimate  the  time  of  chemical  reaction  between  interacting  elements  A  and  B 
Arrhenius  dependence  for  constant  of  the  rate  of  chemical  reaction  may  be  used; 

/nr 

where  A  is  the  frequency  factor.  E  is  the  threshold  energy  of  the  reaction.  E.^,j  shows 

the  decrease  of  the  activation  threshold  due  to  input  of  non-heat  energy,  T  is  the 
temperature. 
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.1  Re-distribution  of  the  released  energy. 

The  energy  of  the  system  (automaton)  released  due  to  the  reaction  is  re-distributed 
(including  its  transformation  into  mechanical  energ},^  of  interaction  with  the  neighboring 

automata)  during  some  certain  time  of  relaxation .  The  value  of  /,-g/^yis  of  the 
same  order  that  the  time  of  transition  to  thermal  equilibrium  state. 

3.  Simulation  of  shock  compression  of  powder  mixture 

To  demonstrate  the  abilities  of  the  MCA  method  to  simulate  shock  compression 
of  powder  mixture  consider  the  Ni— A1  powder  mixture.  Powder  mixture  is  simulated  as  an 
assembly  of  randomly  packed  circle  particles  with  certain  variance  and  appointed  porosity. 
A1  and  Ni  powder  particles  are  randomly  distributed  in  the  mixture.  Each  particle  consists 
of  close  packed  elements  of  the  same  sort.  In  this  task  the  mass  ratio  of  Ni  and  A1  particles 
corresponds  to  stoichiometric  ratio  of  the  reaction  (0.01  and  0.007  respectively).  Size  of  a 
separated  element  was  equal  to  0.001  sm. 


Fig.  2.  Initial  stmcture  of  the  sample. 

The  mixture  is  placed  in  the  hard  box  simulated  by  elements  with  fixed  coordinates 
(Fig.  2).  Special  layer  imitates  plastic  material  of  the  top  of  tlie  mixture.  Shock  loading  is 
simulated  by  setting  up  constant  velocity  to  elements  combining  a  movable  piston.  The 
value  of  loading  velocity  was  amounted  to  1  km/sec.  Further  it  can  be  seen  some  results  of 
simulation  of  shock  loading  of  Ni-Al  powder  mixture.  Fig.  3  shows  evolution  of  the 
structure  of  powder  mixture  under  loading.  One  can  see  fracture  of  Ni  particles  and  filling 
of  the  voids  between  Ni  particles  by  the  Al.  Also  some  reaction  at  the  interface  between  of 
Ni  and  Al  is  seen.  Fig.  4  and  5  show  yelocity  fields  and  trajectories  of  elements  of  the 
simulated  system,  respectively,  at  some  time  steps  during  loading  process. 
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b)  ^ 

Fig.  3.  a)-b)  Structure  of  powder  mixture  at  different  stages  of  loading. 


Fig.  4.  a)-b)  Velocity  fields  of  powder  mixture  at  different  stages  of  loading. 
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Fig.  5.  Trajectories  of  elements  in  some  period  of  calculation 


t=0.0}is 


t=0.25ps 


t=0.5|j.s 


t=0.75)j.s 


t=1.0f.is 


l=1.25ps 


Fig,  6,  Collision  of  aluminum  and  nickel  particles.  Structures  at  different  time  moments 


t=O.O^s 


t=0.5(j.s 


Fig.  7.  Collision  of  aluminum  and  nickel  particles.  Velocities  of  automata. 


Fig.  8.  Collision  of  aluminum  and  nickel  particles.  Velocities  of  automata  at  the  interface 
(t=1.25|as). 

4.  Modeling  of  particle  collision 

Let  consider  the  problem  of  interaction  of  two  particles.  Initially  aluminum  (upper) 
particle  moves  with  velocity  of  800  m/s  and  nickel  one  does  not  move.  The  size  of 


particles  is  .300  pm,  size  of  automata  is  0,3  pm.  Figs,  6  and  7  show  the  results  of 
simulation  at  different  stages  of  the  process. 

In  Fig,  6  one  can  see  the  process  of  mi.xture  of  automata  of  different  particles  at 
the  i'nler-particle  interface.  It  is  ver>'  important  that  MCA  method  allows  to  investigate  the 
process  of  interaction  at  meso-scale  level.  In  Fig.  7  one  can  see  that  the  velocity  of 
aluminum  automata  moving  to  the  crack  in  nickel  particle  is  higher  tlian  automata 
velocities  at  the  interface.  In  Fig.  8  one  can  detail  distribution  of  automata  velocities  at  tlie 
inter-automata  interface.  It  is  important  that  there  is  intensive  shear  deformation  at  tlie 
interface. 

5.  Hot  spot  formation  under  high  rate  shear  deformation 

The  simulation  of  the  Ni-Al  powder  mixture  showed  that  an  intensive 
mechanical' mixing  of  the  components,  caused  by  penetration  effects,  takes  place  under 
shock  compression  of  powder  mixture.  Tliis  leads  to  localized  initiation  of  chemical 
reaction.  To  investigate  the  mechanisms  of  that  mixing  on  atomic  level  special  modeling 
was  carried  out  by  molecular  dynamics  method.  The  simulation  showed  that  a  great 
number  of  “hot  spots”  may  be  foraied  under  considered  conditions.  The  problem  of  hot 
spots  is  very  important  for  studying  the  solid-phase  chemical  reactions.  It  is  necessary  to 
say  that  the  hot  spots  at  certain  critical  size  may  initiate  chemical  reactions  in  the 
material. 

The  hot  spot  is  an  ensemble  of  atoms,  tlie  average  kinetic  energy  of  which 
exceeds  to  a  great  extent  over  the  average  kinetic  energy  of  the  sample.  According  to  the 
Maxwell-Boltzman  distribution  of  velocities  of  atoms,  there  are  always  some  groups  of 
atoms,  tlie  energy  of  wliich  are  higher  than  the  average  one.  However  these  groups  can  not 
be  considered  as  a  hot  spot,  since  they  have  such  an  energy  for  a  short  time  (comparable  to 
a  single  atomic  oscillation).  It  also  should  be  noted,  that  in  the  perfect  crystals  the  hot 
spots  never  appear,  even  though  the  sample  is  subjected  to  heating  or  high  energy  loading. 
Hot  spots  are  not  observed  in  tlie  material  with  structural  defects,  if  it  has  been  relaxed 
and  in  equilibrium. 

The  mechanism  of  formation  of  hot  spots  was  studied  in  the  work  [4]  based  on 
molecular  dynamics  method.  As  a  rule,  the  hot  spots  appear  in  the  regions  witli  high 
density  of  defects  due  to  interaction  with  different  perturbation.  In  the  present  work  tlie 
interaction  of  solitaiy  waves,  appearing  under  shear  loading,  with  vacancy  clusters  was 
investigated.  At  first  the  initial  coordinates  and  velocities  of  atoms  of  a  simulating  system 
are  determined.  Then  the  equations  of  motions  are  calculated  on  the  base  of  known 
interatomic  potentials  and  forces.  Three-dimensional  ciystallite  of  Al,  containing  2 
vacancy  clusters  of  4  vacancies  for  each  was  simulated.  The  calculations  were  carried  out 
using  the  unique  software  “Monster  MD“  [5].  This  softw^are  allows:  1)  to  simulate  the 
behavior  of  materials  with  different  kind  of  structural  defects  under  external  loading;  2)  to 
draw  fields  of  distribution  of  atom  velocities,  stresses,  internal  energies,  local 
temperatures,  radial  distribution  functions  for  different  time  moments,  and  etc.  The 
method  of  molecular  dynamic  is  the  base  of  this  software.  The  structure  of  software  can  be 
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easily  changed  (other  mclhods  of  numerical  integration,  boundar)'  conditions,  c.xternal 
loading,  interatomic  cffcctix  c  potential,  structure  of  sample,  an  other  metal  or  alloy  and  so 
on  can  be  used). 

Three-dimensional  crystallite  of  aluminum,  including  near  9000  atoms.  \\as 
modeled  at  the  temperature  near  to  OK.  Coordinate  axes  were  directed  as:  02  -  along 
<  1 1 1  > ,  OY  -  along  <  2  1  T  > .  OX  -  along  <  0 1  1  >  directions.  Periodic  boundary' 
conditions  were  used  in  OY  and  OZ  directions.  On  the  two  side  surfaces  normal  to  OX 

direction  the  boundary'  condition  were:  K  =  F.  ==  0;  V'.  -  0  ;  -  const,  (wher'c 

is  the  component  of  the  velocities  on  the  left  side  surface  of  the  sample:  V’  is  the 

component  of  the  velocities  on  the  right  side  surface  of  the  sample).  The  interatomic  pair 
interaction  potential  was  calculated  on  the  base  of  pseudopotential  theoiT,  either  as  in 
work  [5]. 

To  initiating  the  solitary'  waves,  for  the  first  time  it  was  discovered  in  [6,7],  a 
shear  loading  was  applied  on  the  sample  for  40  time  steps  (time  step  At-2.42*10'’^  sec). 
Then  propagation  and  interaction  of  solitary  waves  with  the  defects  of  structure  were 
studied. 


Fig.  9.  XY-  projections  of  atomic  velocities  in  perfect  cry'stallite  under  shear  loading  with 
the  rate  of  200  m/sec  at  the  900-th  time  step. 

On  the  first  stage  of  studies  the  mechanical  shear  loading  of  an  ideal  cry'stallite 
was  simulated.  The  results  of  simulation  are  shown  in  Fig.  9.  At  first  the  waves  have  only 
X-component  of  velocities,  they  are  compression-extension  waves  and  propagate  with  the 
speed  close  to  the  sound  velocity  for  this  material.  These  waves  may  be  considered  as  a 
certain  analogue  of  elastic  forerunners.  Then  the  waves  propagating  with  much  less 
speeds  (near  3700m/sec)  are  obser\'ed.  They  are  characterized  by  considerably  greater 
amplitude  and  have  X.  Y  and  Z  components  of  velocities.  The  calculations  showed  that, 
when  changing  the  shearing  rale,  the  propagating  speeds  of  the  solitary'  waves  remain 
constant,  but  the  amplitudes  of  the  waves  herewith  are  changed.  Particularities  of 
propagation  and  interactions  of  such  waves  were  studied  in  [6.7].  It  was  shown  that  they 
may  be  formed  in  heterogeneous  materials  even  at  low  mechanical  loading,  in  particular, 
as  a  result  of  creation  of  “shear  bands”,  microdamages,  at  jump-type  grow  ing  of  cracks 
and  etc. 
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On  the  second  stage  of  studies  the  interaction  of  perturbations,  caused  by  shear 
loading,  witli  vacancy  clusters  were  simulated. 

For  this  purpose  the  shear  loading  with  different  rates  of  200niysec.  50iu/scc  and 
.30in/sec  were  applied  to  the  left  side  surface  of  the  sample.  The  results  of  such  simulation 
are  shown  in  Fig.  10  It  should  be  noted  that  the  hot  spot  arises  in  the  defect  region.  The 
parameters  of  hot  spot  such  as  the  starting  time  of  heating,  the  local  temperature,  the 
duration  of  hot  spot  life  were  discovered  to  depend  extremely  on  the  value  of  shear  loading 
rates.  The  region  containing  vacancy  clusters  will  be  heated  to  higher  temperature,  the 
heating  of  this  region  will  start  at  much  earlier  time  and  duration  of  the  hot  spots  life  will 
be  decreased  by  increasing  the  amplitude  of  solitary  wave  (increasing  of  shear  loading).  It 
follows  to  note  that  the  size  of  hot  spot  will  not  be  changed  under  reducing  the  shear 
loading  from  200m/sec  to  30  m/sec.  This  proves  that  the  nature  of  arising  the  hot  spot  is 
connected  with  elastic  energy  release  which  accumulated  in  the  region  of  the  defect. 


Fig.  10.  XY-  projections  of  atomic  velocities  in  the  sample  with  a  vacancy  cluster  under 
shear  loading  with  the  rate  of  200  m/sec  at  the  time  400-th  steps  . 

The  creation  and  development  of  the  hot  spots  are  accompanied  with  essential 
stnictural  relaxation  in  the  region  of  their  locations.  Moreover  the  most  displacements  are 
observed  towards  OZ  axis  perpendicular  to  the  direction  of  shear  loading  and  the  solitary 
waves  propagation.  During  the  process  of  structural  re-arrangement  the  velocities  of 
separate  atoms  may  reach  very  high  values  (1000  m/sec). 

Thus,  the  energy  of  the  hot  spot  is  defined  by  the  avalanche  process  caused  by 
decreasing  of  potential  energy  under  the  local  structural  re-arrangement  initiated  by 
mechanical  loading.  The  obtained  results  allow,  in  particular,  to  explain  some  effects  of 
mechanical  activation  of  solid-phase  chemical  reaction  of  the  materials  components  [8-10] 
and  to  better  understand  the  response  of  the  materials  under  shear  loading. 

Profile  of  potential  relief  was  shown  to  be  changed.  Conditions  for  low 
temperature  thermo-nuclear  synthesis  were  shown  to  be  reached  in  such  systems. 
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THREE  PEARLS  OF  SHOCK  COMPRESSION  INORGANIC 

CHEMISTRY. 

Yu.Horie 

(North  Carolina  State  University,  Raleigh.  USA) 

Tliere  are  many  problems  in  shock  compression  chemistry  of  inorganic 
materials  that  cry  out  for  solution.  This  paper  is  concerned  with  three  of  the  most 
outstanding  issues:  strong  chemical  initiation,  growth,  self-sustainability. 

A  review  will  be  made  of  recovery  experiments  and  real-time  measurements 
to  infer  or  delineate  controlling  material  responses  that  initiate  strong  chemical 
reactions  at  the  shock  front.  It  will  be  proposed  that  analogous  to  crack  propagation, 
chemical  initiation  is  primarily  controlled  by  local  mechanisms  and  can  not  be 
described  by  global  conditions  such  as  input  energy.  An  approach  based  on  local 
differential  particle  velocity  will  be  discussed  and  an  elementary  initiation  model 
will  be  evaluated  using  such  mixtures  as  Sn/S  and  Al/Ni. 

Chemical  kinetics  at  present  cannot  be  discussed  quantitatively  due  largely 
to  lack  of  measurements.  However,  we  shall  present  model  calculations  of  pressure 
measurements  to  suggest  (1)  that  global  kinetics  has  a  characteristics  time  constant 
on  the  order  of  1  pm  or  less  and  (2)  that  the  amount  of  reaction  at  the  shock  front  is 
limited.  We  surmise  that  this  limit  is  a  result  of  the  fact  that  the  initiation  takes 
place  in  spatially  localized  regions  of  reactant  materials.  Thus,  the  extent  of  fast 
reaction  is  also  controlled  by  the  same  mechanisms  that  initiate  the  reaction. 

Self-sustained  reaction  (condensed  phase  detonation)  is  of  interest  because 
it  provides  unequivocal  proof  of  strong,  ultra-fast  chemical  reaction  at  the  shock 
front,  we  shall  discuss  a  model  calculation  of  self-  sustained  shock-induced  chemical 
reaction  in  a  Ni/Al  mixture,  based  on  the  idea  of  layered  composite  material  to 
induce  and  to  support  the  chemical  reaction. 
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MATHEMATICAL  MODELLING  OF  SHOCK  WAVE  PROCESSES  IN 
CHEMICALLY  REACTING  INORGANIC  POWDERS. 


V.  M.  Fomin.  S. P. Kiselev.  A. P.Alhimov,  O.B. Kovalenko 
(Institute  of  Theoretical  and  Applied  Mechanics  Rus.  Acad,  of  Sci.,  Novosibirsk. 

Russia) 

The  behavior  of  tlie  mixture  of  non-organic  chemically  reacting  powders  in 
shock  waves  is  discussed.  The  influence  of  strength  properties  of  the  materials  and 
inhomogenity  of  deformation  of  the  powders  on  the  process  of  chemical  reactions 
and  formation  amorphous  alloys  at  the  surface  of  the  powder  particles  is  analyzed. 
The  structure  of  the  shock  wave  and  dilatation  wave  is  studied  in  powder  mixture. 
The  method  of  cool  gas-dynamic  spraying  is  considered  for  obtaining  composite 
materials.  The  paper  presents  the  results  of  mathematical  modeling  of  propagation 
of  SHS  stationary  wave  in  Ni-Al  powder  mixture. 


X-DIFFRACTION  INVESTIGATION  CONCERNING  STRUCTURES  OF 
BE,AL,SI,FE,  LIF,S102,KCL  UNDER  DYNAMIC  PRESSURES  FROM 

2GPA  TO  20  GPA. 

V.V.Mokhova.  L. A. Egorov,  A. I. Barenboim,  V.V.Dorokliin 
(Russian  Federal  Nuclear  Center  (VNIIEP),  Arzamas- 16) 

High  velocities  of  deformation  during  effect  of  shock  waves  on  solids  cause  fast¬ 
going  structural  changes  with  generation  of  metastable  intermediate  phases  and 
states.  To  investigate  these  processes  it  is  necessary  to  employ  structural  -  sensitive 
methods  of  solid  physics. 

Currently,  the  only  direct  method  to  study  behaviour  of  solid  crystal  substance 
structures  under  dynamic  compression  is  method  to  record  X-rays  diffraction 
pictures  of  crystal  structures  under  shock  compression. 

The  paper  presents  results  of  X-rays  diffraction  measurements  concerning 
structural  parameters  of  shock  compressed  substances  at  pressures  higher  than 
Hugoniot  elastic  limit  (Be,  Al,  LiF,  Fe+3%Si),  lower  than  Hugoniot  elastic  limit  (Si, 
Si02,  LiF)  and  under  pressure  of  phase  transformation  begining  (KCl). 

Recorded  states  of  shock-compressed  substance  structures  demonstrate  indentity  of 
structural  deformations  at  pressures  higher  and  lower  than  Hugoniot  elastic  limit  as 
well  as  at  pressures  above  the  phase  transformation  point,  which  can  be 
characterized  as  single-axial  defonnations. 
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HIGH-SPEED  CHEMICAL  REACTIONS  AT  SHOCK  LOADING  OF 
CONDENSED  SYSTEMS 

V.A.  Skripnyak.  E  G.  Skripnyak,  E.V,  Shilko* 

Tomsk  Stale  University.  Tomsk.  Russia, 

*  Institute  of  Strength  Physics  and  Materials  Science  SB  Russian  Academy  of  Sciences, 

Tomsk,  Russia 

Introduction 

In  the  case  of  chemical  processes  in  solids  we  deal  with  a  new  class  of 
phenomena,  which  are  different  in  principle  from  the  chemical  processes  in  reacting 
liquids  and  gases.  The  scientific  interest  to  chemical  processes  and  mechanisms  of 
reactions  ongoing  in  solids  is  stimulated  by  3  reasons; 

1.  Possible  difference  in  reaction  products  of  tire  same  initial  components  [  1-3  ]. 

2.  The  realization  of  chemical  reactions  in  the  solid  phases,  which  carmot  be  realized  in 
the  liquid  phases  [1]. 

3.  Difference  in  the  stages  of  chemical  reactions  in  solids  and  liquids.  Solid  phase 
chemical  process  goes  for  a  single  stage,  while  the  reactions  in  liquid  phases  can  proceed 
in  several  stages  with  formation  of  metastable  intermediate  products  [1-3]. 

Solid  phase  chemical  reactions  can  be  used  for  production  of  new  chemical 
products  and  materials  and  to  release  the  energy  with  high  rates  that  lead  to  explosion. 

It  is  well  known,  that  at  usual  conditions  the  chemical  processes  in  solid-phase 
systems  proceed  more  slower,  than  in  liquid  and  gas  phases.  It  is  stipulated  by  slow 
speed  of  diffusion,  ensuring  the  mass  transfer  of  components  in  solid  phase.  At  present, 
several  methods  are  used  to  accelerate  the  chemical  reactions  in  solid  phases  [1-13]. 

The  methods  of  solid-phase  chemical  reactions  acceleration  are  connected  with 
mechanical  loading  of  powder  mixtures.  The  solid  components  are  processed  into  fine 
powders  with  the  grain  size  of  about  10-100  pm.  High  speed  solid-phase  chemical 
reactions  of  decomposition  or  sintering  may  be  realized  in  powder  mixture  under  quasi¬ 
static  mechanical  loading.  This  loading  is  a  combination  of  high  hydrostatic  pressure 
and  intensive  plastic  shearing.  At  laboratory  conditions  high-speed  solid-phase  reactions 
are  obtained  under  loading  in  circuit  rotated  -  Bridgmen  anvil  and  under  mechanical 
extrusions  of  mixture  through  a  head  of  cold  extruder.  The  time  of  the  reactions  is  lO'L 
10"^  sec  [  1-4  ].  The  characteristic  time  of  the  shock-assisted  chemical  reaction  is  of 
about  10'^-  10'^  sec  in  such  powder  mixtures  as  Sn-S,  Mo-S  [  3  ].  Thus,  in  the  both  cases 
the  high  rate  chemical  reactions  take  place.  The  Arrhenious  effective  factor  of 
diffusivity  will  have  the  value  of  about  10  decimal  orders  higher,  than  at  usual 
conditions. 

High  rate  chemical  reaction  can  be  produced  under  liquid-phase  sintering  and 
by  technology  of  self-propagating  high-temperature  synthesis  (SHS)  [5-7]. 

The  time  of  liquid-phase  sintering  of  previously  compacted  powder 
mixtures  of  inorganic  materials  is  of  about  lens  minutes. 
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The  speed  of  combustion  front  propagation  is  4-iO  mm/sec  in  the  mixtures 

(Mo-S , 

Ti-C,  Ti-Si  etc  ).  The  appropriate  characteristic  time  of  the  chemical  transitions  in  SHS 
front 

is  about  1-10  sec.  Direct  observation  oxer  the  SHS  process  in  a  thin  layer  of  mixture 
showed,  that  the  reaction  goes  at  formation  of  thin  layers  of  liquid  phase  [5].  The  rough 
estimation  of  the  particle  heating  time  in  the  reacting  mixtures  gives  10'^  sec.  The 
diameter  of  the  grains  is  less  then  0.5  mm  and  maximum  temperature  corresponds  to  the 
melting  point.  Therefore,  the  reactions  at  SHS  are  not  caused  by  thermal  activation  due 
to  the  heat  transfer  from  the  heated  layers.  It  is  possible  to  solve  this  paradox,  if  to 
assume,  that  the  formation  of  the  melting  zone  is  stipulated  by  the  energy  release  in  the 
zone  during  the  process  of  exothermic  chemical  reaction.  This  process  is  self 
accelerated. 

The  most  effective,  but  not  enough  studied  method  of  acceleration  of  practically 
important  classes  of  solid-phase  reactions  (  decomposition,  restoration  and  synthesis  )  is 
a  shock-wave  loading  of  powder  mixtures. 

The  experimental  study  of  solid-phase  chemical  reactions  under  shock  loading 
are  being  carried  out  for  more  than  30  years  already.  In  Russia  the  basis  of  this  direction 
was  founded  by  G.A.Adadurov,  S.  S.Batsanov,  A.N.Dremin,  Yu.A.Gordopolov, 
A.G.Merzhanov  and  others.  The  overview  of  this  works  was  done  in  [3], 

Up  to  now,  the  studying  of  solid  phase  chemical  processes  under  shock  loading 
are  performed  by  experimental  methods.  At  the  same  time  the  grate  interest  is  exliibited 
to  research  the  solid-phases  chemical  processing  by  mathematical  modeling.  At  present 
there  is  no  united  approach  to  describe  all  kinds  of  high-speed  solid-phases  chemical 
reactions.  Therefore,  it  is  of  great  importance  to  define  both  the  application  fields  of  the 
existing  approaches  to  investigation  of  high-speed  chemical  processes  and  necessity  of 
development  of  new  approaches. 

The  main  question  is  ;  which  of  already  developed  mathematical  models  and 
metliods  are  applicable  to  research  of  high-speed  chemical  processes  and  mechano- 
chemical  phenomena.  In  this  paper  we  shall  discuss  the  mathematical  models  of 
chemical  processing  only  under  shock  loading. 

Models  of  mechano-chemical  processing  under  shock  loading 

Experimental  studies  showed  that  the  phenomena  of  high  rate  solid  phase 
chemical  reactions  are  strongly  differed  from  chemical  processes  in  mixture  of  liquids 
and  gases  [1-14].  The  solid-phase  chemical  processes  are  defined  mostly  by  the 
mechanical  factors  (stress,  intensity  of  plastic  deformation,  uniformity  of  a  stress-strain 
condition  )  and  less  by  thermodynamic  factors  (  such  as  temperature).  That  is  caused  by 
the  difference  between  mechanisms  of  mass  transfer  at  chemical  processes  in  solids  and 
diffusion  mechanisms  in  gases  and  liquids. 

The  basis  for  modeling  mechano-chemical  processes  in  gas  and  liquid  mixtures 
are  the  continuum  mechanics  and  Arrhenious  chemical  kinetics  [  3  ]. 
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As  a  nilc,  a  model  of  mechano-cliejiiical  processes  includes  three  groups  of 
equations  18-10],  The  first  group  consists  of  equations  of  mechanical  conser\'ation  laws 
and  laws  of  phenomenological  thermodynamics.  The  second  group  includes  the 
equations  of  slates  for  multiphase  medium.  Kinetic  equations  of  chemical  processes  are 
united  in  the  third  group. 

The  continuum  models  of  mechano-chemical  processes  in  powder  mixtures 
under  shock  loading  can  include  the  same  groups  of  equations  [8-10],  However,  the 
specific  forms  of  the  equations  may  be  different. 

There  are  some  difficulties  at  the  development  of  the  models  of  solid  state 
mechano-chemical  processing.  Some  of  this  difficulties  are  connected  with  adequate 
describing  of 

-  the  mechanical  behavior  of  porous  powder  mixtures  under  loading; 

-  the  energy  release  at  the  phase  transfomiation  under  loading; 

-  the  change  of  the  solid  phases  concentrations  at  chemical  reactions  under  loading. 

Mechanical  behavior  of  powder  mixture  under  shock  loading  is  described  by  the 
approach  of  multiphase  media  mechanics  [10].  The  solid  phase  chemical  reaction  is 
known  to  be  initiated  under  application  of  hydrostatic  pressure  and  shearing. 

One  of  tile  problems  of  tliis  approach  is  the  necessity  of  construction  of  state  equation  for 
the  mixture,  the  phase  structure  of  which  is  varied  at  chemical  transformation.  It  should 
be  taken  into  account,  that  the  initial  mixture  have  pores,  wliich  are  collapsed  under 
shock  loading. 

Plastic  deformations  are  more  non-uniform  on  meso-level  of  porous  powder 
mixture,  than  in  continuous  solids  under  shock  loading.  Tlierefore  the  chemical 
reactions  are  initiated  in  local  micro-volumes. 

Distribution  of  particle  velocity  in  front  of  shock  waves  was  observed  in  metals 
and  alloys  [  13  ].  Dispersion  of  particle  velocity  behind  the  front  of  the  shock  wave  is  not 
decreased  in  materials  with  essentially  heterogeneous  structure  (for  example  in  steel  and 
high-strengtii  alloys).  This  can  be  explained  by  the  fact,  that  the  shear  stress  relaxation 
keeps  going  behind  the  front  of  tlie  shock  wave,  but  with  the  less  speed.  The  most 
probable  mechanisms,  ensuring  the  shear  stress  relaxation  behind  the  front  is  the 
rotational  mechanisms  of  plastic  flow. 

Metallographical  studies  of  shock-loaded  metals  and  alloys  proved  the  reality  of 
such  mechanisms  [14-15]. 

Two-level  micro-dynamic  model  can  be  offered  for  analysis  of  mechano- 
chemical  phenomena.  This  approach  was  developed  for  studying  the  connection  between 
micro-localization  of  plastic  deformations  and  particle  velocity  dispersion  in  front  of 
plane  shock  waves  in  alloys  and  metals  [16]. 

In  the  framework  of  this  approach  the  mechanical  behavior  of  powder  mixture 
are  considered  on  two  connected  structural  levels.  At  macro-level  the  effective 
mechanical  parameters  are  used.  This  parameters  are  averaged  over  representative 
volume  of  the  mixture. 

The  material  point  of  the  medium  on  meso-scopic  level  is  considered  as  heterogeneous 
media. 
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Sy  stem  of  consen  ation  equations  on  macro-level  includes: 

—  =  -div  u  ,  pu  =diva ,  pE-~(j{pl p)  +  pQ  ,  (1) 

P 

where  p, , .  p  -  initial  and  current  average  mass  density.  ?7  -  particle  velocity,  a  -  stress 
tensor  ,  E  -  specific  internal  energ>%  0  -  rale  of  energy 

The  mechanical  parameter  o*f  state,  mechanical  and  release  in  chemical 
react  ion.  thermodynamic  properties  have  a  distributions  on  meso-level. 

The  initially  two  component  powder  mixture  can  be  considered  as  4  phases 
mixture  after  initiation  of  chemical  reaction.  In  the  simplest  case,  there  are  2  initial 
solid  phases,  1  solid  phase  of  reaction  product  and  1  gas  phase  in  pores. 

The  equation  system  of  phase  mixture  includes 

=  1  ,  P,C-P  ,  j  =  l,4  (2) 

J  =  l,4 

where  Q  -  concentration  of  j  -  phase,  pj  -  mass  density  of  j  phase  in  representative 
volume. 

Density  of  gas  phase  can  be  supposed  to  be  equal  to  zero.  Volume 
concentrations  of  initial  solid  phases  Ci,  C2  and  volume  concentration  of  the  product  of 
chemical  reaction  C3  are  connected  with  stoichiometric  relation 

K,p,  C,+K,p,C,  =K,p,C3  ,  (3) 

where  Kj  -  stoichiometric  coefficients. 

The  specific  internal  energy  can  be  calculated  for  representative  volume 

pE  =  Pj  Ej  . 

(4) 

The  local  themiodynamic  equilibrium  is  assumed  within  representative  volume.  Thus, 

P,E.=-Pj{P,l Pi)+siJ„„,  ;  =  i2 , 

£3  =  -  {p,  /  pP  +  SiJ„„  +  p,  Q,  ,  ,/•  =  3. 

(5) 

where  ^mn  "  components  of  stress  tensor  deviator  and  strain  tensor  in  j-th 

component,  Qj  ~  specific  energy  release  . 

The  internal  energy  of  representative  volume  pE  equal  sum  of  internal  energy  of 
solid  phases  p,  Ej  under  adiabatic  loading  in  front  of  shock  waves. 

The  value  of  can  be  supposed  equal  to  zero  in  hydrodynamic  approach. 

It  is  known,  that  chemical  reactions  of  sintering  arc  initiated  in  powder  mixture  under 
shock  pressure  of  15-  20  GPa  [3,  4]  . 

The  Mie-Gmneisen  equation  of  state  of  solid  phase  can  be  used  [  3,  17  ]; 

P,-P,(P,,,E,  ), 

(6) 
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The  pressures  P,  are  assumed  to  be  equal  to  each  other  in  solid  phases; 

Pi=P2=P..  (7) 

The  models  of  damaged  media  are  popular  for  calculation  of  mechanical 
behavior  of  porous  materials  under  dynamic  loading  [8-10  j 

The  effective  pressure  P  can  be  calculated  on  macro-level  within  the 
framew'ork  of  the  model  of  damaged  media  j  10]. 

P-Pj  (  1-F(C,)) 

(9) 

where  F(C4  )  -  phenomenological  function. 

For  metal  powder  mixtures  it  can  be  approximately  linear  F(C4)  ~  C4/C4 ,  C4  is 
the  constant. 

The  decreasing  of  volume  concentration  of  gas  phase  in  the  pores  of  the  mixture 
are  usually  described  by  evolutionary  or  kinetics  equations  [8-10],  The  main  problem  is 
to  use  the  most  adequate  kinetics  equation  for  the  mixture.  The  simple  kinetic  equation 
can  be  used; 

q=[p-f;(i-f(c,)]/t,p;,  p>  c,>o  , 

C4  =  o,  C4=o, 

(10) 

where  T4  -  characteristic  time  of  the  pores  collapse,  Pq  -  constant  of  mixture. 

It  is  established,  that  the  solid-phase  chemical  transformations  can  be  realized 
both  in  self-consistent  front  of  chemical  reactions  and  in  front  of  shock  waves.  In  the 
first  case  the  speed  of  the  reaction  front  can  be  less  than  the  speed  of  the  shock  wave  [ 
11]  .  Thus,  the  mechanisms  of  chemical  transformation  can  be  different  in  the  dilferent 
solid  phase  mixtures  under  shock  loading.  It  should  be  supposed  that  intensity  and  speed 
of  energy  release  at  solid-phase  chemical  transformations  play  the  important  role  in 
realization  of  this  mechanisms. 

The  energy  release  under  chemical  transformations  ought  to  be  included  in  the 
equation  of  energy  conservation  law  (1).  The  rate  of  specific  energy  release  must  be  also 
allowed  in  (5)  for  chemical  transformation  of  components  at  meso-level. 

The  method  for  calculation  of  specific  energy  release  was  offered  by  S.S. 
Batsanov  for  two-component  mixtures  [3,13,18].  This  method  may  be  recommended  for 
calculation.  If  the  specific  energy  release  Q  is  known,  then  the  rate  of  specific  energy 
release  can  be  calculated  as 

ft  =  2ft  (II) 

The  problem  comes  up,  when  the  mixture  has  three  or  more  reactants.  The  system  of 
equations  (l)-(ll)  is  not  closed.  The  equation  of  chemical  transformation  ( 

Q  =  ,...)  )  ought  to  be  add  to  (l)-(l  1). 

The  change  of  the  phase  concentrations  at  chemical  processing  is  described 
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by  kinetics  equation  of  chemical  reaction.  There  is  no  universal  kinetics  equation  for 
high  rate  solid  phase  chemical  reactions.  Two  questions  have  to  be  solved  for  developing 
kinetic  equations  of  high  rate  chemical  reaction  : 

a) .  What  are  the  conditions  of  loading  which  initiate  high  rate  solid  phase  chemical 
reaction  ? 

b) .  How  the  solid  phase  chemical  transformation  evolves  under  loading  ? 

Initiation  of  high  rate  solid  phase  chemical  reactions 

It  is  necessary'  to  introduce  the  description  of  distribution  of  mechanical 
parameters  into  evolutionary  equation.  It  may  be  constructed  wdth  the  usage  of  the  result 
of  modeling  the  chemical  processes  in  the  mixture  by  molecular  dynamics  method. 
Criteria  of  initiation  and  the  end  of  the  chemical  transformations  can  be  constructed, 
using  parameters  of  phase  deformation  on  meso-level.  The  kinetic  or  evolutionary 
equations  have  to  be  designed  in  terms  of  phase  concentrations  and  mechanical 
parameters  (both  averaged  and  phase  ones). 

The  plastic  deformation  is  the  most  important  factor,  that  controls  of  the 
initiation  and  the  process  of  chemical  reaction.  Collapse  of  the  pores  in  mixture  is 
caused  by  plastic  deformation  of  grains  of  the  mixture.  Thus,  the  effective  bulk  plastic 
deformation  of  macroscopic  medium  is  defined  by  decreasing  volume  concentration  of 
pores  C4.  The  following  approximation  can  be  used; 


where  C  ^  is  the  initial  volume  concentration  of  pores  . 

As  it  was  mentioned  above,  two  kind  of  chemical  processes  are  realized  in 
pow'der  mixture  under  shock  loading.  The  processes  of  the  first  kind  are  shock-initiated 
and  the  second  ones  are  shock-assisted  [3], 

Short  time  of  solid  phase  chemical  reaction  initiation  requires  super  high  speeds 
of  mass  transfer  on  atomic  or  molecular  structural  level.  The  opportunity  to  realize  such 
mechanisms  was  proved  in  [4J.  At  the  same  time,  high-speed  of  solid-phase  chemical 
processes  in  front  of  shock  waves  were  found  to  be  caused  by  mass  transfer  at  the 
distances  of  about  some  pm.  However,  this  mechanism  does  not  explain  chemical 
transformations  in  powder  particles  of  the  mixture  at  the  depth  up  to  100pm.  Thus,  after 
initiation  of  reaction  the  new  mechanisms  of  high-speed  mass  mixing  can  occur. 

Direct  and  indirect  estimation  of  temperature  in  front  of  shock-initiated 
chemical  reactions  show's  that  the  high-speed  solid-phase  chemical  reactions  are  not 
thermally  activated  [11-13].  For  example,  thin  polyethylene  film  in  a  mixture  w'as  not 
melted  after  the  passage  of  the  front  of  exothermic  solid-phase  chemical  reaction  through 
the  mixture  [12].  At  propagation  of  the  shock-initiated  solid-phase  reaction  of 
decomposition  of  bicromat  ammonia  the  introduced  thin  films  of  low-melting  metals 
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(Sn  and  Pb)  remained  not  melted,  while  the  film  made  of  Wood  alloy  was  completely 
melted.  The  speed  of  the  chemical  reaction  front  was  equal  to  1.3  km  /  sec,  that  is  less 
then  the  speed  of  the  shock  wave  propagation  in  bicromat  ammonia. 

The  model  of  solid-phase  chemical  reaction  initiation  has  to  include  the  three  conditions: 

1) .  The  concentration  of  reactant  in  mixture  ought  to  be  more  than  0. 

Cj>0,  .1=1,2,  (13) 

where  Cj  is  the  volume  concentration  of  chemical  components. 

2) .  The  second  condition  can  be  formulated  in  two  forms,  either  for  pressure  or  for 
internal  energy.  The  applied  pressure  has  to  be  higher  than  come  critical  value  P*,  which 
depends  on  the  mixture  porosity,  average  grain  size  and  distribution  of  grain  size. 

P>P’.  (14) 

Specific  internal  energy  of  each  reactant  Ej  ought  to  be  liigher  than  critical  value  Ej*. 

E,>e;  ,j=1,2.  (15) 

The  specific  internal  energy  of  components  Ej  are  different  under  adiabatic  loading  in 
the  shock  waves. 

3) .  The  third  condition  is  connected  with  the  specific  mechanisms  of  high  rate  solid 
phase  chemical  reactions.  The  effective  plastic  deformation  of  the  medium  ought  to  be 
increased. 


where 


dsf  >0, 

>  0, 


(16) 


(dsL)^3l 


is  the  increment  of  intensity 


_1  P 

of  effective  plastic  deformation,  Q  s^k  "  increment  of  the  bulk  effective  plastic 
deformation  (12). 

The  energy  release  under  exothermic  chemical  reaction  leads  to  the  rise  of  non¬ 
uniformity  of  temperature  and  stress  fields  and  initiates  the  rotational  modes  of 
deformation  on  different  scale  levels.  Within  this  framework  the  two  following 
experimental  phenomena  can  be  explained;  1)  the  dependence  of  critical  of  pressure  of 
initiation  of  solid-phase  reaction  on  the  size  of  powder  particles  and  2)  the  decrease  of 
the  critical  pressure  at  introduction  of  inert  powder  inclusions  into  tlie  mixture  [2], 


Kinetic  equations  of  chemical  reactions 

The  rate  of  solid-phase  reactions  under  shock-wave  loading  is  not  constant.  Last 
year  experiments  carried  out  in  the  Institute  of  Chemical  Physics  (Moscow)  showed,  that 
the  chemical  reactions  can  be  initiated  in  the  very  front  also.  Shock  front  structure  was 
studied  in  reacting  two-component  mixtures  Sn-S,  Mg-S,  Ni-Al  by  pyrometric  method 
[3].  The  time  of  solid-phase  reactions  is  about  50-200  nano-sec  for  these  mixtures  [3, 
18]. 

Application  of  Arrhenious  kinetic  equation  in  computer  modeling  of  solid- 
phase  chemical  processes  needs  to  define  the  effective  kinetic  coefficients  in; 
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(17) 


^=F({C,},{k}). 

d  t 

where  t-  ihe  time,  F  -  a  fmictjon.  C,  -  the  specific  volume  of  component  1  ,  k  -  the 
parameter  connected  with  the  speed  of  chemical  reaction. 

Specific  volume  of  reaction  product  C3  can  to  calculate  by  (2)-(3),  (17).  The 
kinetic  equation  (17)  can  be  written  in  linear  fornt  in  the  simple  case: 


d  t 


The  value  of  k  is  not  constant  in  the  phase  mixture  on  meso-level.  It  is  the  complex 
characteristics  of  the  loading  process  and  tlie  average  size  of  powder  particles; 

lc=k{E,U,T},  (19) 

where  E  -  an  activation  energ>%  U  -  the  work  of  stresses,  T  -  the  temperature  in  the 
energy'  units. 

The  equation  (19)  can  be  used  in  the  form: 

k  =  V(,exp[-(E+U)/T]  ,  (20) 

where  Vq  k  lO’  ^  c'’  . 

The  characteristic  time  of  exothermic  synthesis  reaction  in  solid-phase  systems 
Sn-S  ,  A1  -Fe203  etc.,  measured  by  pyrometric  measurements  of  intensity  of  radiation  on 
the  border  of  the  reacting  mixture,  is  of  about  0. 1-0.5  10'®  sec  [  3  ].  Thus,  the  rale  of 
solid  phase  chemical  reaction  in  powder  mixture  under  shock  loading  can  be  of  3-8 
decimal  orders  higher  than  tlie  rate  of  similar  reaction  in  a  liquid  phase.  The  estimation 
based  on  the  speed  of  front  propagation  of  shock-initiated  solid-phase  chemical  reaction 
gave  the  same  values  [11]. 

The  calculation  results  of  plane  shock-wave  loading  of  stoichiometric  Ti-Al 
powder  mixture  with  application  (l)-(U)  ,  (12)  -(16)  and  (18)  are  shown  in  Fig.l  and 
Fig.2  .  Profiles  of  the  shock  wave  front  are  shown  corresponding  to  the  time  moment  0.4 
p,sec  after  the  impact. 

The  initiation  of  chemical  reaction  is  not  realized  in  front  of  the  shock  wave  in 
Fig.l.  The  condition  of  formation  of  two-wave  structure  of  shock  front  is  defined  by  the 
kinetics  of  pore  collapse  (9)-(10).  If  the  chemical  reaction  is  initiated  in  front  of  the 
shock  wave,  the  front  struclure  is  transformed.  Shock -assisted  chemical  reaction  leads  to 
acceleration  of  shock  wave  front  propagation  and  some  rise  of  the  pressure.  Results  of 
calculation  of  shock  wave  structure  in  that  case  are  shown  in  Fig.2.  Non-uniformity  of 
distributions  of  internal  energies  and  strains  of  solid  phases  are  risen  on  meso-level  of 
mixture  after  beginning  of  chemical  reaction.  The  chemical  conversion  is  not  full  in 
stoichiometric  mixture  (  A1  ,Ti  ->  Ti  AI  )  after  shock  wave  loading  with  amplitudes  of 
13  GPa.  The  specific  volume  of  TiAl  slowly  grows  with  the  distance  from  the  impact 
surface.  The  chemical  reaction  is  finished  when  specific  volume  of  pores  comes  to  zero. 
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Conclusion 


The  multi-level  model  of  solids  can  be  used  for  modeling  chemical  processing 
of  solids  under  shock  loading.  Multi-phase  mcchamcai  models  of  reacting  gases  or 
liquid  mi.xtures  are  rough-casts  for  the  solids.  The  distributions  of  mechanical 
parameters  on  mcso-level  of  solid  mixture  includes  the  information  about  possibilib'  of 
the  initiation  and  the  process  of  chemical  reactions. 
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Fig.l.  Structure  of  the  shock  wave  front  in  stoiheometric  Ti-Al  mixture  with  initial 
specific  volume  of  pores  0.1.  Time  after  impact  corresponds  to  0.4  10'*^  sec.  The 
chemical  reaction  was  not  initiated.  I  -  pressure  P,  2  -specific  volume  of  pores  (C4). 
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Fig.2.  Structure  of  the  shock  wave  front  in  stoiheometric  Ti-Al  mixture  with  initial 
specific  volume  of  pores  equal  to  0.1.  Time  after  impact  corresponds  to  0.4  10'^ 
sec.  The  chemical  reaction  was  not  initiated.  1  -  pressure  P  ,  2  -  specific  volume 
of  pores  (C4),  3  -  specific  volume  of  TiAl  (C3),  which  are  created  in  the  shock 
wave  front ,  4  -  point  of  chemical  reaction  initiation. 
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INTERFEROMETRIC  MEASUREMENT  OF  THE  PARTICLE  VELOCITY 
DISPERSION  IN  DYNAMICALLY  LOADED  SOLIDS 

Yu.l.  Mescher>’akov 

Institute  of  the  Mechanical  Engineering  Problem 
Russian  Academy  of -Sciences.  Saint-Petersburg 

1 .  Introduction. 

Mechanics  of  deformed  solid  is  traditionally  based  on  the  concept  of  mean 
stresses  and  mean  strains.  However,  the  processes  of  strain  localization  cannot  be 
understood  and  described  without  introducing  an  intermediate  scale  level  of 
deformation  -  the  so-called  mesoscopical  level  [1].  To  date  there  are  several  scientific 
schools  in  Russia  which  developed  and  advocate  mesomechanical  approach  to  the 
strain  localization  problems  [2,3]. 

As  for  dynamic  plasticity,  the  only  possibility  to  watch  a  behavior  of 
elementary  carriers  of  deformation  (ECD)  at  the  mesolevel  is  to  introduce  into 
consideration  particle  velocity  distribution  function  (VDF).  In  practice,  however,  its 
determination  seems  to  be  very  difficult  if  possible  at  all.  Therefore  the  measuring  of 
statistical  moments  of  VDF  becomes  of  great  importance.  In  the  shock-wave 
experiments  as  a  rule  only  the  first  statistical  moment  of  VDF  is  measured,  namely  the 
time  profile  of  average  particle  velocity.  In  case  of  uniaxial  strain  loading  this  profile  is 
registered,  for  example,  with  capacity  or  interferometric  gauges.  In  [4,5]  the  technique 
for  determination  of  the  second  statistical  moment  of  VDF,  particle  velocity  dispersion 
<  V  (t)>  ,  has  recently  been  developed.  For  different  materials  it  has  experimentally 
been  shown  that  particle  velocity  dispersion  does  not  remain  to  be  constant  along  tlie 
average  velocity  profile  V  (t)  but  currently  changes.  For  example,  for  the  steady  plastic 
fronts  maximum  value  of  dispersion  occurs  at  the  middle  of  front.  As  for  the  unsteady 
plastic  fronts  the  behavior  of  dispersion  is  not  so  definite.  It  may  often  increase  along 
the  plastic  front  up  to  pulse  plateau. 

Microstructure  analysis  of  specimens  after  shock  loading  shows  that  there 
exists  a  certain  correlation  between  particle  velocity  dispersion  at  the  mesolevel  and 
strain  localization  processes.  It  was  found  that  kind  of  kinematical  mechanism  at  the 
mesolevel  (translational,  i.e.  shear  banding,  or  rotational)  depends  on  the  value  of  tlie 
velocity  dispersion  during  the  shock  wave  passage. 

The  objective  of  this  paper  is  to  show  that  kind  of  kinematical  mechanism 
depends  not  only  on  the  longitudinal  component  of  velocity  dispersion  but  on  the 
transverse  component  of  dispersion  as  well. 

2.  Mesolevel  and  statistical  approach. 

Technique  for  determination  of  mean  particle  velocity  in  a  shock  loaded 
material  (the  first  statistical  moment  of  VDF)  is  based  on  the  interferometric  measuring 
Doppler  frequency  shift  of  laser  radiation  when  the  latter  drops  on  the  free  surface  of 
shock  loaded  target.  In  case  of  the  so-called  velocity  interferometer  mean  particle 
velocirf  V  (t)  relates  to  fringe  number  N(t)  with  a  well-known  relation.  \6\: 

=  (1) 

2t 
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where  x  is  the  interferometer  delay  time  When  the  particle  velocities  at  the  free  surface 
of  target  are  equal  to  each  other,  so  that  the  scattering  of  particle  velocities  are  absent, 
interference  signal  amplitude  is  constant  and  fringe  frequency  changes  in  accordance 
with  changing  the  free  surface  velocity.  If  however,  within  laser  beam  spot  at  the  free 
surface  of  target  there  exists  a  scattering  of  the  particle  velocity,  amplitude  of 
interference  signal  decreases  due  to  widening  of  laser  radiation  spectrum  after 
reflecting  from  the  target  [4,5], 

Until  now,  when  speaking  about  particle  velocity  dispersion  under  conditions 
of  dynamic  deformation  we  did  not  give  concrete  expression  for  their  kind  and  space 
scale.  It  is  not  hard  to  estimate  their  typical  scale  by  using  the  principle  of  interference. 
Measuring  the  particle  velocity  dispersion  is  based  on  concept  that  different  particles  at 
the  free  surface  of  target  within  laser  beam  spot  give  different  Doppler  frequency  shifts. 
Effective  interaction  of  separate  particle  and  laser  radiation  is  possible  only  if 
dimensions  of  particles  are  more  greater  as  compared  to  laser  radiation  wave  length  Xo 
(approximately  ten  times).  On  the  other  hand,  uniform  widening  of  radiation  spectrum 
is  possible  only  in  case  if  within  laser  beam  spot  Dsp  at  the  free  surface  of  target  there  is 
sufficient  quantity  of  similar  particles  (at  least  about  hundred).  From  here  the  following 
estimate  for  dimension  of  particle  can  be  accepted: 

X.O  dp  Dgp 

As  a  rule,  being  focused  on  the  free  surface  laser  beam  has  a  transverse  dimension 
about  100  pm  and  laser  radiation  length  equals  «  pm  from  where  the  value  of  particle 
approximately  equals  dp  =  10pm.  That  dimension  in  accordance  with  generally 
accepted  in  physics  of  plasticity  classification  belongs  to  the  so-called  mesoscopical 
scale  level. 

Note  that  the  expression  which  relates  interference  contrast  I  and  particle 
velocity  dispersion  <  V  (t)>^ ,  [5] 


I  =  exp 


(3) 


It  has  been  obtained  under  assumption  that  VDF  is  equilibrium,  i.e.  has  a  symmetrical 
maxwellian  form.  This  is  true  only  for  the  steady  plastic  fronts.  In  general  case  of 
unsteady  fronts  in  [8-  9]  alternative  technique  for  determination  of  the  particle  velocity 
distribution  width,  has  been  developed.  It  is  grounded  on  the  phenomenon  of  the  so- 
called  "fluctuative  decay”  of  waves  in  a  structure  nonuniform  media.  Consider  this 
approach  in  details. 

In  accordance  with  classification  suggested  by  G.E.  Duvall  [7],  there  exist 
three  kinds  of  wave  decay:  i).  geometrical  decay,  ii)  hydrodynamical  decay,  iii) 
maxwellian  decay.  The  first  kind  of  decay  relates  to  space  geometry  of  shock  wave.  For 
the  spherical  wave  decay  is  proportional  to  1/  r^  .  for  the  cylindrical  wave  decay  is 
proportional  to  1/r  and  for  the  plane  wave  decay  does  not  depends  on  propagation 
distance. 
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Hydrodynainical  kind  of  deca>  relates  to  the  circumstance  that  back  front  of 
compressive  pulse  propagates  with  the  higher  ^elocily  than  the  first  front  since  the 
former  moves  in  a  compressed  material.  For  the  short  three-angle  pulses  this  results  in 
"cutting"  upper  part  of  pulse.  Lastly  the  so-called  "ma.xwellian"  decay  is  stipulated  by 
dissipative  processes  in  a  medium. 

Beside  the  enumerated  kinds  of  \\a\e  decaw  which  can  c.xist  in  uniform 


medium,  in  nonuniform  medium  an  additional  kind  of  wave  decay  appears,  which 
results  from  the  nonuniform  character  of  the  particles  motion  -  the  so-called 
"fluctuative  "  decay.  Its  nature  can  be  understood  from  the  statistical  consideration  of 
the  wave  propagation  process  in  a  structure-nonuniform  media. 

Aforepresented  estimate  shows  that  space  scale  of  structure  nonunifonnity 
attributes  to  mesoscopical  level.  Investigations  in  physics  of  plasticity  which  were 
carried  out  during  the  last  ten  year  have  shown  that  one  of  the  most  typical  features  of 
mesolevel  is  the  space  and  charge  heterogenization  of  dislocation  structure  on  the 
distances  of  the  order  of  free  run  of  dislocations.  Simplified  scheme  of  heterogenization 
of  dislocation  ensemble  is  presented  in  Fig.  1 .  It  is  seen  that  majority  of  dislocations 
from  the  mesovolumes  V  is  subdivided  in  sign  and  gathered  into  narrow  layers  AV  in 
vicinity  of  the  mesovolume  boundaries.  These  ranges  have  abundant  dislocation  density 
Ap.This  means  that  these  layers  can  be  considered  as  superdislocatins  having  some 
summary  Burgers  vector  B.  Ensemble  of  similar  particles  are  known  to  have  a  strong 
long-range  interaction,  which  in  turn  superimposes  the  specific  features  on  their 
statistical  behavior.  In  other  words,  mesovolume  with  polarized  dislocation  structure 
can  be  represented  as  a  quasiparticle  with  effective  charge  and  therefore  has  a  more 
higher  radius  of  interaction  as  compared  to  initial  quasineutral  state  in  which  Ap  and 
accordingly  Aq  are  small.  This  permits  for  description  of  mesoscopical  plasticity  to  use 
the  statistical  approach  based  on  application  of  Fokker-Plank  equation  for  mesoparticle 
velocity  distribution  function: 


df 

- + 

at 


V 


dx 


m  av 


f)  +  f--^(D,  f) 

)  2  gy  2  V  2  ! 


(4) 


Here  m  is  the  mesoparticle  mass,  F  is  the  external  force. 

Diffusion  coefficients  Di  and  D2  characterize  change  of  VDF  due  to  random 
interactions  of  mesoparticles  with  each  other.  The  first  diffusion  coefficient 

a(Av) 

D,  =  — ^ -  ,  or  dynamic  friction  coefficient,  characterizes  the  change  of  mean  flow 

velocity  due  to  mutual  interaction  of  mesoparticles.  Here  <Av>  is  averaged  on  the 
particle  ensemble  change  of  mean  velocity  owing  to  interaction  of  particles  with  each 

e(Av) 

Other.  Thus,  the  value  D,  =  — ^ — -  has  a  sense  of  deceleration  and  being  multiplied 

by  mesoparticle  mass  defines  the  value  of  fluctualive  friction  force  directed  opposite  to 
the  driving  force.  As  a  matter  of  fact,  this  value  defines  the  fourth  kind  of  decay  which 
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Cfin  be  classified  as  a  "fluctuatiNc  decay",  li  disappears  when  parlicai  velocity 
dispersion  becomes  negligible. 

The  second  diffusion  coefficient  in  the  Fokker-PIank  equation 

e(AV  AV) 

^2  = - - defines  a^  eragcd  on  (he  mesoparticles  ensemble  tlie  rate  of  change 

of  particle  velocity  dispersion.  As  it  has  been  showm  in  (8|.  in  case  of  long-range 
mutual  interaction  of  particles  between  diffusion  coefficients  there  exists  the  following 
relation: 

A=/r+^(A)=i>r'+p^:  (5) 

where  is  the  invariable  part  of  friction  coefficient  related  to  polarization  of 
particle  ensemble.  It  follows  from  equation  (5)  that  fluctuative  decay  depends  on  the 
rate  of  changing  the  second  diffusion  coefficient.  This  conclusion  is  of  great 
importance  for  understanding  the  nature  of  fluctuative  decay.  In  a  sense  it  means  that 
change  of  the  mesoparticle  velocity  dispersion  ,  for  example  due  to  nucleation  of  shear 
bands,  results  in  fluctuative  braking  of  mesoparticles,  i.e.  decreasing  the  pulse 
amplitude. 

In  Fig.  2  two  different  free  surface  average  velocity  profiles  obtained  for  two 
kinds  of  aluminum  alloys  are  presented  to  compare  different  features  of  material 
depending  on  the  dispersion  behavior.  Both  profiles  has  been  obtained  under  the  same 
impactor  velocity  of  160  m/s.  Under  symmetrical  shock  loading  when  impactor  and 
target  have  identical  acoustic  impedance  there  is  following  relation  between  impactor 
velocity,  particle  velocity  and  free  surface  velocity:  V^p  =  2Vp  =  . 

It  is  clearly  seen  from  Fig.2  that  while  in  aluminum  AMg-6  maximum  pulse  amplitude 
corresponds  to  impactor  velocity  in  aluminum  alloy  V-95  free  surface  velocity  at  the 
plateau  of  pulse  achieves  only  110  m/s,  so  that  decrease  of  the  mean  particle  velocity 
equals  approximately  to  50  m/s.  This  decreasing  occurred  owing  to  of  catastrophic 
nucleation  of  shear  bands  under  free  surface  velocity  of  110  m/s  .  Thus,  abrupt 
increasing  a  mesoparticle  dispersion  results  in  fluctuative  decay  of  average  velocity 
pulse. 

Note,  further,  that  at  the  back  front  of  compressive  pulse  a  number  of 
interference  fringes  is  again  restored.  This  means  that  fluctuative  decay  at  that  part  of 
pulse  is  absent.  From  the  point  of  view  of  microstructure  kinetics  this  means  that 
random  mutual  motion  of  mesoparticles  relative  each  other  is  stopped  since  a 
development  of  shear  bands  results  in  relaxation  of  internal  microstresses.  Thus, 
fluctuative  decay  proves  to  be  reversible  process.  At  the  first  pulse  front  the  energy 
related  to  motion  of  mesoparticle  with  mean  velocity  pumps  into  energy  of  their  chaotic 
motion.  During  the  pulse  plateau  relaxation  of  internal  stresses  leads  to  decreasing 
dispersion.  This,  in  turn,  results  in  disappearing  the  fluctuative  decay.  Schematically 
the  process  of  changing  a  mean  velocity-  can  be  represented  in  form  of  changing  the 
shape  of  particle  velocity'  distribution  function  as  it  is  shown  in  Fig.3.  Under  free 
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surface  velocit>’  of  1 H)  m/s  a  desymmelrizalion  of  Ihe  VDF  happens  wliereas  to  the  end 
of  pulse  plateau  VDF  again  restores  to  s>'mmclrical  shape. 

|l  should  be  noticed  that  total  rc\’ersibilily  of  transition  from  s>  mmetrical  VDF 
to  nonsymmetrical  one  and  vice  versa  is  possible  onl>  for  the  elastic  interactions  of 
mcsoparticles.  In  the  case  of  iinclastic  interaction  the  energy  loss  at  the  e.xpense  of 
unelastic  dissipative  processes  takes  place.  For  c.xamplc,  de\  elopmcnt  of  shear  bands  is 
exclusively  dissipative  process  resulting  in  local  heating  of  material.  It  should  be 
remembered,  however,  that  all  what  related  to  dissipative  losses  is  already  taken  into 
account  in  a  "max^vellian"  decay.  Fluctuative  decay  is  the  pumping  of  energy  from  the 
mesoparticle  motion  with  mean  velocity  into  chaotic  motion  and  vice  versa,  while 
"maxwellian"  decay  is  non-reversible  energy'  scattering.  In  a  certain  sense,  fluctuative 
decay  can  be  considered  as  a  latent  energy  box  from  which  under  specific  conditions 
the  energy  of  chaotic  motion  suddenly  transforms  into  energy  of  directed  motion  of 
particles,  which  can  often  result  in  fracture  of  material.  The  same  reason  can  lie  in  a 
unexpected  onset  of  earthquake. 


3.  Kinematics  of  straining  at  the  mesolevel. 

In  this  connection  the  question  arises  how  the  particle  velocity  distribution 
function  and/or  its  statistical  moments  are  related  to  the  kind  of  kinematical 
mechanism  of  dynamic  straining.  One  possibility  to  analyze  this  coupling  seems  to  use 
the  continual  theory  of  dislocation  keeping  in  mind  that  instead  of  dislocation  it  is 
implied  the  mesoparticle  having  a  summary  Burgers  vector  equaled  to  dislocation 
charge  of  mesoparticle.  Complete  equation  system  can  written  as  follows,  [9j ; 


M 


dV,  da. 


dx,  - 


di  dX, 

, 

■  a 


Im 


(6) 

(7) 


Here  ^  is  the  mesoparticle  velocity,  Wik  are  components  of  elastic  distortion, 

{A 


U  is  the  total  displacement  (elastic  plus  plastic)  and  Jj,  are  the  components  of 
mesoparticle  flow  tensor.  Kinematics  of  plastic  straining  can  be  inferred  from  the 
continuity  equation  (7).  By  using  this  equation  it  may  be  written; 


dX, 


d^ 

a 


■  + 


{■4-4) 


(8) 


Since  we  are  interested  only  in  kinematics  of  plastic  straining  elastic  rotations  will  be 
neglected.  Then  the  rate  of  plastic  rotation  is  determined  by  difference  of 
antisymmetrical  components  of  flow  tensor: 


CO  = 


a-',  3\ 


Uy  dx,. 


=  4-  -  4 


(9) 
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In  accordance  wiih  dcfinilioii  of  ^o^\  tensor  \  ia  particle  \  eIocitA  distribution  function 
|10|: 

~  ^ tJ (1^)) 

Its  components-for  the  coiicrcte  rotation  along  coordinate  X^i  are: 

■^21  =  j(f'i/2i  -VJs^dV-,  (11) 

■^2  -^2.  =jy>{/n  +/„Vn-  +f,,)dv 

Here  f  are  the  components  of  tensor  mesoparticle  velocity  distribution  function.  [10]  ; 

7n  /,2  /u' 

/=  /:,  /22  fn  ;  fii  are  the  screw  components 

V31  fsi  f33,y 

4  arc  the  edge  components. 

It  is  seen  from  Eq.(ll)  that;  i)  rotational  motion  of  medium  in  plane  is 

provided  by  motion  of  screw  components  of  VDF  along  axis  X^ ,  ii)  rotational  motion 
achieves  maximum  value  when  motion  of  edge  components  at  plane  X^  X^  is 
suppressed,  iii)  when  the  densities  of  screw  and  edge  components  of  VDF  are 
identical,  difference  between  antisymmetrical  components  of  flow  tensor  equal  to  zero; 

‘^12  “  '^21 

which  means  that  rotational  motion  of  medium  is  absent  and  there  is  only 
translational  motion. 

Let  K3  =  Vjj  is  the  component  of  mesoparticle  velocity  along  the  wave  propagation 
and  =  V2  =  the  components  in  transverse  direction.  If  the  densities  of 

screw  and  edge  dislocations  arc  identical  one  obtains; 

=  (12) 

Let 

Suppose  then  that  longitudinal  and  transverse  mean  mesoparticle  velocities  are  equal 
to  each  other.  Here  A  V^j  and  A  are  the  total  mesoparticle  velocity  decrease  in 

longitudinal  and  transverse  direction  due  to  fluctuative  decay  mechanism.  They 
relate  to  the  first  diffusion  coefficient  D  as  follows: 

M\  =  \lX'di  A(',=J /;,(//  (13) 

B>  using  Eq.  (5)  one  obtains: 
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iy‘ 

Ar,=-~- 


(14) 


where  D:  is  tlie  second  diffusion  coefficicnl  of  the  Fokkcr-Plank  equation.  Thus,  tlic 
rotation  rate  can  be  c.xprcsscd  in  -the  following  form: 


ZJ 


=j/T) 


jyi  ix 

V  ^  //  ^  J  J 


cJV  (15) 


It  is  seen  that  intensity  of  rotation  motion  depends  on  the  difference  between 
variation  coefficient  in  longitudinal  and  transverse  directions.  For  the  constant 
mesoparticle  density  this  equation  can  be-written  in  the  form: 

zu-pB{s^^  (16) 

where 


AffAF 

V~ 

is  the  variation  coefficient. 

It  can  be  concluded  experimental  technique  to  be  used  must  provide  the 
simultaneous  measuring  of  both  longitudinal  and  transverse  dispersion  at  the  single 
action  of  dynamic  straining. 
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Electrical  response  of  (Sn+S)  reactive  heterogeneous 
system  on  shock  wave  effect 

S.S.  Nabatov 

(InstiliUc  of  Chemical  Physics  (Chcrnogolovka).  Rus.  Acad  of  Sci..  Chernogolovka. 

Russia) 

Introduction 

Investigations  of  electrophysical  phenomena  in  inorganic  materials 
understanding  of  the  mechanism  of  shock  induced  transformations  in  heterogeneous 
systems.  Electrical  methods  attract  attention  and  have  promising  perspectives  because 
they  can  be  used  in  dynamic  experiments,  restricted  by  microsecond  time  interval  [1], 
as  well  as  in  reeovery  experiments  with  an  unlimited  registration  time  [2,3], 

In  the  present  report  a  new  method  of  investigation  of  the  (A+B->AB)  type 
reactions  at  shock  compression  is  described.  In  this  case  A  is  metal  component,  B  - 
dielectric  component,  and  AB  semiconductor  compound.  The  method  idea,  first 
formulated  in  [4],  is  based  on  the  well  known  fact  that  emf  generated  at  shock 
compression  of  metals  (1-10  mV)  is  much  lower  than  that  for  semiconductors  (100- 
1000  mV).the  experimental  realization  of  the  idea  consists  in  the  registration  of  the 
electrical  response  of  shock  compressed  mixture  to  the  effect  of  shock  have  reflected 
from  electrode  of  the  experimental  device.  It  is  supposed  that  the  shock  compressed 
mixture  will  generate  sufficiently  large  signal  under  the  effect  of  reflected  stock  wave  if 
semiconductor  compound  is  appeared  behind  the  front  of  the  first  shock  wave. 

The  method  has  been  proved  with  the  system  Sn+S->SnS  in  the  pressure 
range  10-30  GPa  (first  compression).  It  has  been  found  out  that  the  electrical  signals 
registered  at  the  mixture  compression  by  the  reflected  wave  correspond  to 
semiconductor  state  of  the  sample  when  the  pressure  of  tlie  first  shock  wave  exceeds 
15Gpa.  Preliminary  conclusion  has  been  made  that  the  signals  can  be  explained  be  the 
formation  of  the  semiconductor  compound  SnS  behind  shock  wave  front. 


Experiment  and  results 


Shock  wave  experiments  were  carried  out  with  the  use  of  monolithic  samples 
made  from  elemental  powder  mixture  corresponding  to  the  stoichiometiy  of  the  SnS 
compound.  The  size  of  tin  particles  was  about  100  micrometers.  The  scheme  of  tlie 
experimental  arrangement  is  illustrated  in  Fig.  1.  The  sample  with  a  sickness  of  1-2  mm 
and  10  mm  in  diameter  is  sandwiched  between  the  aluminum  driver  plate  and  the  thick 
copper  electrode.  The  sample  and  electrode  are  placed  inside  teflon  ring.  The  mi.xture 
IS  subjected  to  a  plane  shock  wave  compression  with  the  help  of  tlie  high  explosive 
charge  (60  mm  in  diameter)  which  is  in  contact  with  driver  plate. 
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The  pressure  of  the  first  shock  wave  PI  and  the  pressure  of  the  wave  reflected  from  the 
electrode  P2  w'cre  calculated  on  the  base  of  the  standard  shock  impedance  matching 
method  using  the  known  pressure  in  the  dri>  er  plate  and  the  Hugoniots  of  Cu.  A1  and 
the  system  Sn+S.  The  Hugoniot  of  the  system  Sn^S  is  |5] 

D=2. 13+I.62*U.  km/s:  po=4.74  g/cc.  where  D  is- shock  ^•eloclty.  U  -  particle  \'elocity 
and  po  -  initial  density. 

Tlie  voltage  V  of  the  electrical  signal  was  recorded  by  a  fast  digital  oscilloscope  with 
the  input  resistance  Re=75  Ohm.  Since  tlie  sample  resistance  is  small  compared  to  Re 
the  voltage  V  is  equal  to  the  emf  generated  by  the  shock  wave  when  it  propagates 
betw'een  the  electrodes. 

As  it  has  been  reported  m  [6]  SnS  sample  generates  at  shock  compression 
signals  of  characteristic  voltage  200-300  mV.  The  response  form  is  shown  qualitatively 
in  Fig.2.  where  tl  -  time  of  shock  wave  entrance  into  a  sample  t2  -  time  of  shock  wave 
arrival  at  the  second  electrode.  The  system  Sn+S  does  not  reveal  semiconductor 
properties  at  PI  about  10  GPa.  However  at  more  higher  pressures  the  reactive  mixture 
really  behaves  as  a  semiconductor  substance.  This  fact  is  illustrated  by  experimental 
oscillograms  presented  in  Fig. 3.  The  upper  record  (Fig.  3a)  was  obtained  at  P  1=2 1.0 
GPa.  The  first  shock  wave  generates  weak  positive  signal  on  the  level  about  20  mV.  At 
the  moment  of  the  first  wave  reflection  from  the  electrode  (t2)  the  negative  impulse 
with  the  amplitude  more  than  250  mV  arises.  During  the  reflected  wave  propagation 
through  the  mixture  the  response  changes  its  polarity.  Oscillogram  obtained  at  P  1=25.5 
GPa  and  P2=39.6  GPa  (Fig.  3b)  is  qualitatively  similar  to  previous  one  with  the 
exeption  of  tlie  short  positive  impulse  at  the  beginning  of  the  record.  It  is  important  to 
note  that  this  impulse  arises  before  the  shock  wave  reflection  from  the  electrode. 

Discussion 

Within  the  scope  of  the  qualitative  model  of  the  method  it  is  impossible  to 
evaluate  the  extent  of  the  chemical  transformation.  But  there  are  indirect  evidences 
showing  that  the  mixture  does  not  react  completely  behind  the  shock  front.  The  above 
mentioned  short  impulse  at  the  beginning  of  the  oscillogram  in  Fig.  3b  is  an  example 
of  such  evidence.  The  rise  of  the  impulse  before  reflection. of  the  first  wave  from  the 
electrode  can  be  explained  be  the  following  way.  Due  to  partial  chemical  interaction  the 
contacts  between  tin  particles  were  destroyed  and  the  mixture  acquired  semiconductor 
properties  behind  the  front  of  the  first  shock  wave!  Thus  the  compressed  part  of  tire 
mixture  can  respond  as  a  semiconductor  element  on  a  thermal  effect.  The  thermal 
affect  can  be  caused  by  a  full  completion  of  the  reaction  in  some  localized  region  («hot 
spot»)  of  the  compressed  mixture.  To  our  mind  there  is  analogy  between  the  sharp 
signal  in  Fig.  3b  and  signals  that  have  been  observed  during  registration  of  the  reaction 
Sn+S->  SnS  in  recovery'  experiments  [3j.  The  e.xpcrimental  arrangement  for  electrical 
measurements  in  recoveiy-  experiments  is  illustrated  schematically  in  Fig.  4,  The 
sample  under  investigation  is  in 
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contact  with  tlie  semiconductor  element  generating  llicrmo  emf  signal  due  to 
temperature  difference  between  its  sides.  I’he  thermo  emf  signal  obtained  m  the 
e.xpcriment  witli  the  system  Sn+S  is  shown  in  Fig  5.  The  pressure  in  this  expenment 
was  16  GPa  (multiple  compression).  One  can  sec  short  impulses  caused  b>  hot  spots 
bounding  with  semiconductor  (SnS)  element  surface  But  this  case  chemical  reaction 
takes  place  in  millisecond  time  interval. 

It  is  believed  that  the  method  developed  will  find  w  ide  applications  m  the  field 
of  shock  wave  chemistry. 
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IMPACT  INDUCED  SOLID  STATE  METAL/METAL  REACTIONS=^ 


Diana  L.  Woody  and  Jeffen  J  Davis 
Naval  Air  Warfare  Center  Weapons  Division 
Chine  Lake.  California  93555-6001 

ABSTRACT 

This  paper  discusses  e.xpcriinental  results  from  an  effort  conducted  to  discern 
the  basic  mechanism  of  reactions  in  porous  metal/metal  compositions  under  rapid  plastic 
flow  conditions.  Small-scale  impact  tests  w^ere  performed  on  various  intermetallic 
mixtures:  3CuO  +  2Al  Fe203  +  2A1.  Ni  +  Af  and  5Ti  +  3Si.  The  addition  of 
polytetrafluoroethylene  (Teflon)  to  the  metal/metal  mixture  has  been  demonstrated  to 
affect  the  extent  of  the  reactions.  Real-time  emissiviw  of  the  reacting  materials  was  used 
to  discern  the  exothermic  reactions  occurring  under  rapid  flow'  conditions. 

INTRODUCTION 

Metal/metal  and  metal/metal  oxide  materials  have  exhibited  the  ability  to 
generate  highly  exothermic  reactions  capable  of  occurring  in  the  same  time  frame  as  , 

detonations  in  conventional  energetic  materials.’’^  Until  recently,  these  materials  w'ere  | 

assumed  to  initiate  exclusively  under  shock  conditions.  How'ever,  work  done  by  Woody 
and  Davis^'’’^  with  materials  under  mechanical  impact  conditions  as  well  as  by 
Nestereenko  and  Meyers  wfth  materials  under  explosive  compaction  conditions,  has 
demonstrated  the  role  that  rapid  plastic  flow  can  play  in  the  reactions  of  these  mixtures. 
Horie  has  modeled  their  shock  reaction  based  upon  a  critical  plastic  flow.^ 

Small-scale  experiments  have  been  performed  to  observe  the  effect  of  rapid 
plastic  flow  produced  by  impact  reactions  on  porous  metal/metal  compositions.'  This 
paper  explores  the  effect  of  the  addition  of  polytetrafluoroethylene  (Teflon)  to  several 
metal/metal  and  inetal/metal  oxide  mixtures  including  3CuO  +  2A1,  Fe203  +  2AI,  Ni  + 

Al,  and  5Ti  +  3Si.  Reactions  such  as  those  found  in  high-temperature,  self-propagation 
synthesis  (SHS)  have  been  observed  under  impact  conditions.  Shear  due  to  impact- 
induced  plastic  flow  plays  an  important  role  in  the  rapid  energy  release  of  these  types  of 
materials.  The  plastic  flow'  induced  reaction  process  is  comparable  to  shock  assisted 
reactions  and  SHS.  A  two-color  infrared  detector  w'as  used  for  real-time  emission 
measurements  of  the  reacting  materials, 


•  Approved  for  public  release;  distribution  is  unlimited. 

EXPERIMENT 

A  drop  weight  impact  machine  was  used  to  induce  a  plastic  flow  in  the  test 
samples.  The  impact  machine  used  for  these  experiments  consisted  of  an  anvil  an 
accelerated  guided  drop  w'eighl,  a  base,  and  a  release  triggering  device.  The  impact 
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machine  is  described  fijlJy  in  another  publication  ^  Elastic  shock  cords  were  used  to 
accelerate  the  drop  weight  to  obtain  impact  A^elocities  of  13  m/s.  At  these  velocities,  the 
plastic  floe  has  been  measured  to  be  approximately  60  m/s.  The  impact  of  the  drop 
weight  on  the  anvil  as  planar  to  AAithin  2  mrad.  The  samples  used  for  the  tests  weighed 
0.2  g  and  \Aere  in  loose  powder  from  prior  to  impact.  Because  the  samples  were  loose 
powders,  their  initial  porosity  was  not  measured  The  light  emanating  from  the  impacted 
sample  came  from  its  surface. 

Teflon  was  mixed  with  the  metal  mixtures  described  above  by  the  mortar  and 
pestle  method.  The  specifics  of  the  materials  are  listed  in  Table  1. 


TABLE  I.  Characteristics  of  Materials  Used  in  This  Study. 


Material 

Manufacturer 

Material  Data 

Titanium  (Ti) 

CERAC,  Inc. 

325  mesh,  99.5%  purity,  20pm  or  less 

Silicon  (Si) 

CERAC,  Inc. 

325  mesh,  99.5%  purity,  10pm  or  less 

Aluminum  (Ai) 

CERAC,  Inc. 

325  mesh,  99.97%  purity,  5  - 15  pm 

Teflon 

DuPont 

7 A,  35  pm,  60%  crystallinity 

Iron  Oxide  (Fe203) 

CERAC,  Inc. 

325  mesh,  99.97%  purity,  >15  pm 

Copper  Oxide 
(CuO) 

;  CERAC,  Inc. 

325  mesh,  99.97%  purity 

Nickel  (Ni) 

1  CERAC,  Inc. 

i  325  mesh,  99.5%  purity,  10  -  20  pm 

Garnet  Paper 

Norton,  Inc. 

:  180A  Garnet  A51 1 

The  two-color  infrared  detector  consisted  of  a  HgCdTe  element  juxtaposed  to  an 
InSb  element.  Each  element  measured  0.101  by  0.101  cm  with  am  active  area  of  0.010 
cm^.  The  elements  were  housed  in  a  liquid-nitrogen-cooled  Dewar  and  kept  at  an 
operating  temperature  of  77  K.  The  InSb  element’s  spectral  response  was  from  2  to  5 
pm.  The  HgCdTe  element  was  capable  of  detecting  wavelengths  from  5  to  12  pm.  The 
signal  from  each  infrared  detector  element  was  transmitted  as  a  voltage  tlirough  an 
initial  voltage  amplifier  and  the  transferred  to  a  LeCroy  digital  oscilloscope. 


RESULTS  AND  DISCUSSION 

The  reaction  of  tlie  materials  to  impact  was  quantitatively  defined  from  tire 
infrared  detector’s  signals.  The  extent  of  the  reaction  was  quantified  from  the  area  under 
the  infrared  emission  curve.  The  most  extensive  reactions  in  these  experiments  were 
defined  as  the  largest  peak  emissions  and  the  shortest  time-to-peak  emission  registered 
by  the  two-color  infrared  detector.  The  peak  emission  and  lime-to-peak  emission  for 
some  of  the  experiments  are  given  in  Table  II, 
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TABLE  II.  Results  Obtained  with  the  Infrared  Detectors 


for  Some  of  the  Compositions  Tested. 


Material  and 

Composition 

Radiatio 
n  From 
InSb 
Detector, 
mV 

Radiation 

From 

HgCdTe 

Detector. 

mV 

Time-to- 

Peak 

Emission,  ps 

5Ti  +  3Si 

29 

1.41 

2950 

(99%)  5Ti  +  3Si 
(1%)  Tefton  (35  urn) 

16 

6.4 

2950 

(90%)  5Ti  +  3  Si 
(10%)  Teflon  (35  urn) 

3547 

55 

56 

(80%)  5Ti  +  3  Si 
(20%)  Teflon  (35  pm) 

3664 

68.9 

56 

(80%)  Ti 

(20%)  Teflon  (35  pm) 

2890 

64.5 

31.5 

3CuO  +  2A1 

3250 

164.1 

16.25 

(90%)  3CuO  +  2A1 
(10%)  Teflon 

3090  ' 

163.8 

27 

(75%)  3CuO  +  2A1 
(25%)  Teflon 

3000 

158.1 

44 

(50%)  3CuO  +  2A1 
(50%)  Teflon 

62.5 

10.2 

750 

Fe.Os  +  2A1 

22.7 

3.9 

2300 

(99%)  FezOg  +  4A1 
(1%)  Teflon 

3280 

94.5 

3075 

(90%)  Fe203  +  4A1 
(10%)  Teflon 

3500 

154.7 

36 

Ni  +  A1 

34.8 

10.9 

1340 

(90%)  Ni  +A1 
(10%)  Teflon 

63.3 

7.7 

810 

(50%)  Ni  +A1 
(50%)  Teflon 

62.5 

10.2 

750 

Note:  Percentages  are  by  mass. 


Qualitative  signs  of  reactions  consisted  of  such  parameters  as  a  visible  light 
emission,  an  audible  signal,  and  characteristics  of  the  recovered  simple’s  surface.  The 
most  exothermic  samples  exhibited  a  \'isible  flash  and  sustained  burning  after  impact  up 
to  1  to  3  seconds  in  duration:  charring  was  also  found  on  the  reco\'ered  sample.  It  was 
obsen^ed  that  the  burning  duration  could  be  varied  b)'  changing  the  percentages  of 
Teflon  added  to  the  mixtures.  The  results  obtained  with  the  infrared  detectors  are  given 
in  Table  II. 
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5Ti  +  3  Si 


As  shown  in  Table  II.  impact  of  the  5Ti  -i-  .iSi  mixture  generated  a  small 
emission,  as  recorded  by  the  infrared  detectors.  The  small  emission  coupled  with  the 
longer  time-to-pcak  emission  indicated  a  heating  of  (he  impacted  metals  due  to  dynamic 
compaction  rather  than  a  chemical  reaction.  X-ray  diffraction  of  the  recovered  samples 
reinforced  our  conclusion  that  reaction  did  not  take  place  in  that  system.  The  addition  of 
Teflon  to  5Ti  +  3 Si  considerably  increased  the  peak  emissions  recorded  by  the  infrared 
detectors  and  reduced  the  tine-to-peak  emission  b)'  tow  orders  of  magnitude  (from 
milliseconds  to  microseconds).  A  detailed  study  has  been  performed  on  the  effect  that 
varying  the  percentage  by  weight  of  Teflon  added  to  various  5Ti  +  3Si  mixture  has  upon 
the  exothermic  reactions  observed.^ 

3C.uO  +  2A1 

An  extensive  exothermic  reaction  was  observed  for  the  neat  3CuO  +  2A1 
mi.xture.  The  reaction  was  characterized  by  a  large  infrared  emission  from  the  two 
detector  elements,  a  visible  flash,  and  sustained  burning.  The  spectrometer  registered 
strong  emission  lines  for  the  neat  mixture.  The  3CuO  +  2 A1  mixture  was  the  only 
material  in  our  study  that  reaction  without  die  addition  of  Teflon.  The  addition  of 
polymers  to  the  metal/metal  mixtures  has  been  demonstrated  to  affect  the  extent  of  the 
reactions.  As  can  be  seen  in  Table  II,  increasing  the  percentage  by  weight  of  Teflon 
added  to  the  3CuO  +  2A1  mixture  progressively  decreased  the  peak  emission  and 
increased  the  reaction  time  registered  by  the  infrared  detectors.  It  was  also  shown  that 
varying  the  percentages  of  the  materials  in  the  composites  could  determine  whether  a 
material  is  more  likely  to  produce  a  sustained  burning  reaction.  Studies  were  also 
performed  in  which  the  mixture  containing  90%  (3CuO  +  2A1)  and  10%  Teflon  was 
pressed  into  pellets  to  measure  the  effect  of  porosity  on  reaction  initiation.  The  porosity 
of  the  material  was  a  significant  factor  in  determining  the  extent  of  reaction  and 
exothermic  output  of  these  metal/metal/polymeric  mixtures.  The  current  study 
demonstrated  that,  when  the  mixture  of  90%  (3CuO  +  2A1)  +  10%  Teflon  was  pressed 
into  pellets,  the  %  TMD  had  a  significant  effect  on  the  reactivity  of  tlie  material.  The 
observable  exothennic  sensitivity  decreased  when  the  metal/metal/polymeric  mixtures 
was  pressed  into  a  pellet.  These  data  are  detailed  in  another  publication.^*^ 

Fc^O,  +  4A1 

As  show  in  Table  II,  the  Fe203  +  4A1  mixture  generated  a  relatively  small 
emission  upon  impact.  These  results  were  similar  to  those  obtained  upon  impact  of  the 
5Ti  +  3  Si  mixture.  The  addition  of  Teflon  to  the  Fe203  +  4A1  mixture  considerably 
increased  the  peak  emissions  upon  impact.  Increasing  the  percentages  of  Teflon  added  to 
the  Fe203  +  4 A1  mixture  gad  a  significant  effect  on  the  extent  of  the  exothermic  reaction 
and  the  signature  of  the  infrared  emission  curve.  The  tiine-to-peak  emission  was 
decreased  two  orders  of  magnitude  (from  milliseconds  to  microseconds)  for  the  mixture 
containing  10%  Teflon. 
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Ni  +  Al 


Table  II  shows  that,  although  the  addition  of  Teflon  to  the  Ni  +  A1  mixture 
increased  the  infrared  emission  recorded  by  the  detectors  and  decreased  the  reaction 
time,  the  overall  increase  in  exotlrermic  release  was  not  as  substantial  as  that  obsei^  ed 
during  the  impact  of  the  5Ti  3 Si/Teflon  and  Fe.O.^  +  4Al/Tenon  mixtures.  The 
qualitative  signs  of  the  reaction  mentioned  in  the  preceding  paragraphs  of  this  paper 
were  not  present  these  impact  conditions. 

CONCLUSIONS 

The  results  of  these  experiments  show  that  the  variation  in  the  amount  of  Teflon 
added  to  the  metal/metal  and  metal/metal  oxide  mixtures  can  significantly  alter  the 
extent,  initial  time  of  reaction,  and  the  duration  of  the  exothermic  reaction  upon  rapid 
shear  conditions.  The  addition  of  Teflon  also  influences  the  likelihood  of  a  material  to 
produce  a  sustained  burning  reaction.  The  5Ti  +  3  Si,  and  Fe203  +  4A1  mixture  exhibited 
relatively  similar  responses  to  the  addition  of  varying  percentages  of  Teflon.  Under  the 
impact  conditions  discussed  in  this  paper,  the  highly  exothermic  reaction  of  2  A1  +  3CuO 
did  occur  under  plastic  flow  conditions  without  the  addition  of  Teflon.  However,  when 
Teflon  was  added  to  the  3CuO  +  2A1  mixture,  the  exothermic  reactions  decreased  in 
intensity.  Overall,  the  most  exothermic  reactions  were  observed  for  the  most  porous 
mixtures.  This  gives  an  indication  that  shear  due  to  rapid  plastic  deformation  plays  an 
important  role  in  the  initiation  of  the  metal/polymeric  mixtures  under  impact  conditions. 
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Overview  ;  BOUNDARY  INSTABILITIES  AND  PHASE 
TRANSFORMATIONS  SOLID  MEDIA  UNDER  SHOCK  LOADING 

A.L.Mikhailov,  O.B.  Dremiov.  V.N.  German,  A. I.  Lebedev,  V.A.  Raevski 
(Russian  Federal  Nuclear  Center  (VNIIEP),  Arzamas- 16,  Russia) 


The  phenomena  of  instabilities  t3T3ical  of  boundaries  dividing 
heterogeneous  (stratified)  liquids  or  gases  under  conditions  of  high-velocity  flows 
are  well-known:  Kelvin -Helmgoltz  instability  (slip  instability),  Rayleigh-Taylor 
(gravitational)  instability,  Rjchtmyer-Meshkov  instability  (i.e.  instability  of 
Rayleigh-Taylor  type  under  pulse  effect).  Phenomena  of  such  a  kind  are  less 
studied,  if  we  deal  with  solid  media,  for  example,  metals,  when  strength  and 
viscous-elastic  characteristics  of  media  have  detennining  influence  on  instability. 

The  paper  presents  review  of  fulfilled  experiments  to  investigate 
developments  of  three  instability  t^qies  mentioned  above  on  boundaries  dividing 
heterogeneous  metals  under  effect  of  shock  waves  with  high  intensities  (up  to  50- 
100  GPa). 

The  paper  analyses  relation  between  instability  development  and  the 
strength  as  well  as  phase  transformations  in  media  under  investigation.  Tlie  paper 
illustrates  some  methods  to  stabilize  growth  of  disturbances.  Experimental 
methods  and  diagnostics  are  described. 

Attention  is  played  in  particular  to  investigation  of  phase  transformations 
under  shock  loads  observed  in  the  similar  experimental  set-up. 


MECHANISTIC  PROCESS  INFLUENCING  SOLID  -STATE  SHOCK 
CHEMISTRY  AND  SHOCK  SYNTHESIS  OF  MATERIALS. 

N.N.Thadhani 

(Scool  of  materials  Science  and  Engineering.  Georgia  Institute  of  Technology,  Atlanta. 

Georgia  USA) 

Shock  chemistry  and  its  implementation  in  synthesis  of  materials  requires  the 
development  of  a  fundamental  understanding  of  the  complex  mechanistic  processes  occurring 
during  shock  compression  of  powders.  This  understanding  can  be  obtained  by  combined 
analysis  of  time-resolved  measurements  and  characterization  of  microslructure  of  materials 
obtained  from  shock  recovery  experiments.  In  this  presentation,  we  will  describe  our  recent 
research  work  perlbrmed  using  nanosecond-resolution  time-  resolved  pressure  measurements 
and  detailed  characteri2;ation  of  the  microstructure  of  recovered  materials  in  several 
intermetallic-forming  (e.g.,  Ti-Si  and  Ti-Al )  and  ceramic-forming  (Ti-C  and  Ti-B  ) ,  highly 
reactive,  powder  mixtures. 

Results  of  experiments  performed  at  and  below  the  tlireshold  conditions  for 
occumence  of  reactions  in  the  systems  will  be  shown  that  the  results  of  these  experiments 
clearly  demonstrate  the  unique  effects  of  the  shock -compression  response  of  powders,  which 
are  in  turn  influenced  by  not  only  the  shock  loading  conditions,  but  also  the  morphological 
characteristics  (size  and  shape)  of  individual  powder  particles,  the  volumetric  distribution  of 
the  reactant  powders  in  a  mixture,  and  the  dissimilarities  in  the  mechanical  properties  of  the 
reactants.  For  example  in  the  case  of  Ti-si  system,  powders  of  10-25  pm  particle  size  are  able 
to  plastically  deform  and  flow  (as  revealed  by  post-mortem  raicrostructural  characterization), 
resulting  in  an  intimately  mixed,  highly-activated,  dense-packed  configuration,  which  readily 
undergoes  chemical  reaction  in  the  fiigh-  pressure  shock  state.  Evidence  of  shock-  induced 
chemical  reaction  in  these  powders  is  obtained  by  increases  in  wave  sped  measured  using 
time-  resolved  pressure  measurements.  In  contrast,  powders  of  coarse  particle  sizes  undergo 
fracture,  and  the  fragments  particles  remain  entrapped  without  mixing  with  each  other,  thereby 
remaining  in  an  unreacted  state.  Powders  of  very  fine  particle  sizes  also  show  limited  mixing 
and  inability  to  mix  and  react  due  to  effects  of  agglomeration  of  particles,  as  observed  by  post¬ 
mortem  microstmctural  characterization,  and  as  evidenced  by  time-  resolved  measurements. 
Observations  of  similar  effects  on  reaction  propensity  will  also  be  presented  for  the  Ti-Al,  Ti- 
B,  and  Ti-C  powder  mixtures.  In  general,  it  will  be  shown  that  in  these  various  powder 
mixtures,  the  propensity  for  reaction  can  be  explained  on  the  basis  of  the  ability  of  individual 
reactant  particles  to  undergo  deformation  and  flow,  or  fracture  and  dispersion.  Thus,  conditions 
which  inhibit  the  deformation  (or  Fracture)  and  subsequent  mixing  of  reactants  exhibit  a 
reduced  propensity  for  reaction,  irrespective  of  the  heat  of  reaction  of  the  system. 

For  applications  relevant  to  synthesis  of  materials,  it  will  be  shown  that  the 
mechanochemical  characteristics  also  provide  a  means  of  control  for  “  shock -induced” 
reactions  (  occurring  in  the  shock  state  )  or  '‘shock-assisted”  reactions  (  occurring  during  post- 
thermal  treatments  due  to  shock-  activation  effects).  Thus,  products  with  tailored  stoichiometry 
and  microstructure  can  be  synthesized,  using  either”  shock-assisted”  processes. 
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INVESTIGATION  OF  PHASE  FORMATION  UNDER  THE  EFFECT  OF 
SHOCK  WAVE  IN  THE  V-HF  SYSTEM. 


A. V. Kolesnikov 

(Institute  of  Chemical  Physics(  Chernogolovka),  Rus.  Acad. of 
Sci.,Chernogolovka,  Russia) 


‘Phase  diagram  of  the  V-Hf  system  has  been  investigated  in  detail  at 
equilibrium  conditions.  The  existence  only  of  V2Hf  intermetalide  as  well  as 
solid  solutions  on  the  base  of  initial  components  have  been  established. 

The  use  of  explosions  as  a  preparative  method  permits  to  hope  for 
receiving  of  new  compounds,  which  can  be  of  interest  in  the  chemistry  of 
metallic  superconductors  and  metastable  states. 

Explosion  experiments  have  been  performed  with  the  use  of  cylindrical 
recovery  capsules,  subsequent  X-ray  phase  analysis  of  the  recovered  specimens 
being  carried  out. 

As  the  result  of  the  investigations  two  new  compounds  of  V2Hf3  and 
VHf3  compositions  has  been  found  out  in  the  V-Hf  system  besides  V2Hf 
known  before.  The  compound  V2Hf3  has  sufficiently  wide  region  of 
homogeneity  (50-70  At.  %  H^. 

V2Hf3  rentgenographical  data  permit  to  consider  its  structure  as  the 
disordered  variant  of  the  intermetallic  structure  of  Fd3m  space  group.  The 
disordered  structure  has  been  identified  on  the  base  of  cubic  sell  with  the 
parameters;  a=4,591  Angstrem,  z=l.  X-ray  density  10,5  g/cm3,  Fm3m  space 
group. 

The  compound  VHf3  has  been  identified  on  the  base  of  hexagonal 
primitive  sell  with  the  parameters;  a=4,999  Angstrem,  c=8,151  Angstrem,  X-ray 
density  11,04  g/cm3,  z=2,  Pb3mmc  space  group. 

Thus  the  use  of  the  techniques  of  shock  wave  synthesis  leads  to  an 
essential  extension  of  preparative  methods.  In  particular,  in  the  case  of  the  V-Hf 
system  it  resulted  in  the  synthesis  of  two  new  metastable  compounds.  The 
necesary  condition  of  shock  wave  synthesis  is  the  formation  of  melt. 


CHEMICAL  REACTIONS  IN  CONDITIONS  OF  INTENSIVE, 
CONTROLLED  HIGH-STRAIN  RATE  SHEAR  DEFORMATION 


M,A.  Meyers  ,V.F.  Nesterenko,  ,L  LaSalvia  and  H.B.  Chen 
(  University  of  California,  San  Diego  Xa  Jolla,  USA). 

The  controlled  high-strain-rate  shear  bands,  having  thicknesses  of  5-20 
microns,  were  generated  in  heterogeneous  reactive  porous  materials  (  Nb+Si,  Ti+Si 
)  using  axially  symmetric  experimental  configurations,  based  on  the  radial  collapse 
of  a  thick-walled  cylinder  (  ‘Thick-Walled  Cylinder^’  method).  The  experimental 
method  generated  overall  strains  up  to  100  and  strain  rates  of  approximately  = 
10000000  sec-1  inside  shear  bands.  Changes  in  particle  morphology,  melting,  and 
regions  of  partial  reaction  were  observed  only  inside  shear  bands  for  Nb+Si  system. 
At  the  same  conditions  of  high-strain  rate  deformation,  for  Ti-Si  system  the  reaction 
being  initiated  inside  shear  bands  propagated  tlirough  the  whole  sample.  The 
threshold  conditions  were  established  for  global  chemical  reaction  in  the  volume  of 
Ti-Si  mixture. 


SHOCK  INDUCED  CHEMICAL  REACTIONS  IN  TI-C  AND  TI-B 
POWDER  MIXTURES 

YrM-Fyodorov,  Yu.AJvchenko,  Yu.A.Gordoi)olov 
(Institute  of  Structural  Macrokinetics,  Rus.  Acad.of  Sci.,  Chemogolovka,  Russia  ) 

The  experiments  of  shock-wave  loading  of  exothemiic  mixtures  Ti-C  and 
Ti-B  of  stoichiometric  composition  prove  that  at  relatively  low  pressures  (about 
IGPa)  initiation  of  the  chemical  reaction  occurs  with  complete  transfonnation  of 
initial  rnixmires  into  products  TiC  and  TiB2.  Increase  of  pressure  or  elongation  of 
time  results  first  in  cessation  of  the  chemical  reaction  initiation.  However,  at  10- 
20GPa  partial  transformation  of  initial  components  into  products  TiC  and  TiB2 
takes  place,  and  at  higher  pressures  complete  transformation  is  observed.  Such 
behavior  of  the  mixtures  under  study  can  be  explained  by  the  change  of  mechanisms 
of  chemical  transformation  due  to  the  change  of  shock-wave  conditions. 

The  experiments  prove  that  chemical  transformation  in  the  system  Ti-C  at 
pressures  above  1 1  GPa  occurs  in  the  shock  wave  front  (at  least,  partially). 

The  results  obtained  can  be  used  in  further  investigations  of  mechanisms 
and  kinetics  of  chemical  transfonnations  in  Ti-C  and  Ti-B  systems  upon  shock- 
wave  loading. 
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special  report 

SOFTWARE  BASED  ON  MOVABLE  CELLULAR  AUTOMATA  METHOD 
ORIENTED  TO  SIMULATE  RESPONSE  OF  MATERIALS  AND  STRUCTURES 
UNDER  HIGH  RATE  LOADING 


S.G.Psakhie,  A.  Yu.Smolin,  A.I.Dmitriev, 

E.  KShilko,  S.  Yu. Korostelev,  S.  V.  Alekseev 

Institute  of  Strength  Physics  and  Materials  Science,  Russian  Academy  of  Science,  Tomsk, 

Russia. 


1.  Main  features  of  the  software 

The  MCA-Compaction  software  is  basically  developed  for  simulation  of 
powder  compacting  processes  including  chemical  reactions  under  the  contract  of  Sandia 
National  Laboratories,  USA,  Contract  No  AM  2957B,  The  mam  code  of  MCA- 
Compaction  is  written  in  "C"  programming  language  for  OS  UNIX  System  V  and  is  an 
implementation  of  MCA  Method  [1].  The  software  also  includes  the  utility  MCA- 
Analyzer  to  overview  the  obtained  results  in  on-line  regime  as  well  as  during  the  run  of 
MCA-Compaction.  MCA-Analyzer  is  written  in  "C++"  programming  language  using 
the  XI 1  graphic  interface.  The  MCA-Compaction  software  was  tested  on  SUN  and 
DEC  workstations  and  also  in  CRAY  Y-MP  supercomputer. 

To  run  MCA-Compaction  the  user  have  to  prepare  some  necessary  data. 
Namely  one  have  to  describe  the  materials  characteristics,  the  sample  and  the  parameters 
of  the  task.  Materials  characteristics  describe  mechanical,  physical  properties  and 
parameters  of  phase  and  chemical  transfonnations  and  are  stored  in  materials_db.dat 
file.  The  information  about  the  samples,  including  chemical  compounds,  size  and 
geometry'  of  the  sample,  the  size  of  the  automata  etc.,  is  specified  in  *.smp  filcs.^The 
parameters  of  interaction  among  automata  depending  on  their  types  is  specified  in  .pai 
file.  And  such  parameters  as  the  time  of  calculations,  time  step  of  integration,  type  of  the 
task  (total,  chemical  or  mechanical  calculations  only),  loading  conditions  (velocity  or 
force  value  on  boundaries  and  initial  ones)  etc.,  are  specified  in  *.tsk  files. 

l\/jQ/^_Analyz6r  is  developed  for  X1 1  Window  System  and  utilizes  standard 
XI  ib  library.  This  program  has  main  window  with  the  menu  and  allows  to  display  in 
separate  windows  the  following  characteristics; 

•  structure  of  the  simulated  sample  (automata  are  presented  as  circles), 

•  rotation  angles  of  the  automata. 

•  temperatures  of  the  automata, 

•  velocity  field  of  the  automata. 

•  trajectories  of  automata  between  two  any  time  moments. 

•  net  stnicture  of  the  interacted  pairs  of  the  automata. 

•  scale  for  displayed  results  (in  this  window  user  may  change  the  scales). 

Every'  separate  yvindow  has  the  buttons  for  changing  the  time  step  and  storing  the  results 
in  GIF  format. 


74 


Beside  it  user  has  an  opportunity  to  set  a  number  of  "gages"  in  the  sample,  in  this 
case  MCA-Com paction  will  write  specified  parameters  of  "gage"  automata  to  text  file 
task.dat.  And  then  user  can  \  iew  the  plots  of  result  parameters  versus  time  or  other 
parameter  by  standard  visualization  software. 


2.  Demonstration  of  Possibilities 

To  demonstrate  the  main  possibilities  of  MCA-Com paction  software  the  result 
of  simulations  of  two  problems  are  shown.  The  first  one  is  tension  of  a  sample.  The  size  of 
the  sample  is  0.29x0.53  mm.  The  loading  of  constant  rate  of  50  m/s  is  applied  to  the  top 
of  the  sample.'  In  Fig.  1  the  structures  of  the  sample  arc  shown  for  different  time  moments 
In  Fig.  l,a  the  initial  structure  is  presented,  in  Fig.  l,b  one  can  see  the 


a  be  d  e  f  g 

Fig.  1.  Results  of  modeling  of  sample  fracture  under  tension  a,bx  -  the  case  of  ordinary' 
material;  d.e.fg  -  the  case  where  automata  of  the  sample  are  allowed  to  rotate  without 
breaking  tlieir  bond. 
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beginning  of  fracture  and  in  Fig.  l.c  the  sample  is  broken.  In  the  frame  of  the  MCA 
Method  we  can  varv'  parameters  of  the  interaction  among  automata  and  allow  them  to 
rotate  relatively  each  other  without  breaking  linkage  (bond).  The  results  of  simulation  of 
tension  the  sample  consisting  of  automata  with  these  properties  are  presented  in 
Fig.  l,d,eXg. 

The  second  demonstration  is  bending  of  the  sample  under  tangential  pulse.  The 
size  of  the  sample  is  the  same  as  previous  one.  The  initial  velocity  of  pulse  applied  to  the 
top  of  the  shmpte  is  1 70  m/s.  The  velocities  of  automata  of  the  sample  at  different  time 
moments  are  presented  in  Fig.  2.  One  can  see  that  plasticity  of  the  material  cause 
dissipation  of  initial  pulse  and  at  the  end  the  shape  of  the  sample  is  changed. 


a  bed 

Fig.  2.  Results  of  modeling  of  the  sample  bending  under  tangential  pulse  loading. 


3.  Shock  Chemistry 

At  first  we  consider  the  problem  of  interaction  of  two  particles.  Initially 
aluminum  (upper)  particle  moves  with  velocity  of  800  m/s  and  nickel  one  does  not  move. 
The  size  of  particles  is  300  pm,  size  of  automata  is  0.3  pm.  Figs.  3  and  4  show  the  results 
of  simulation  at  different  stages  of  the  process. 

In  Fig.  3,d,e,f  one  can  see  the  process  of  mixture  of  automata  of  different 
particles  at  the  inter-particle  interface.  It  is  very  important  that  MCA-Compaction 
allows  to  investigate  the  process  of  interaction  at  meso-scale  level.  In  Fig.  4  one  can  see 
that  the  velocity  of  aluminum  automata  moving  to  the  crack  in  nickel  particle  is  higher 
than  automata  velocities  at  the  interface. 

The  following  example  demonstrates  modeling  of  shock  compaction  of  powder 
mixture.  The  results  of  modeling  are  shown  in  Figs.  5  and  6.  In  the  stmeture  windows 
smaller  particles  are  aluminum  ones  and  bigger  -  nickel  particles,  light  gray  automata 
compose  the  flyer. 
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Fig.  3.  Collision  of  aluminum  and  nickel  particles.  Structures  at  different  time  moments. 


Fig.  4.  Collision  of  aluniinum  and  nickel  particles,  Velocities  of  automata 
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Fig.  5.  Results  of  modeling  of  shock  compaction  of  powder  mixture,  a  -  the  structure 
window  of  MCA-Analyzer  with  initial  structure;  b,c  -  computed  stnicture  at  t=25  ps  and 
t=35  ps  respectively. 


b 

Fig.  6.  Results  of  modeling  of  shock  compaction  of  powder  mixture,  a  -  field  of  velocities 
(t=5  ps);  b  -  field  of  temperatures  (1=35  ps). 
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4.  Fracture  Processes  Under  Compression 


-  Sintered  powder  materials  (PM) 

Consider  the  sintering  PM  material  based  on  AI-Cu  alloy  under  compression. 
Special  experiments  was  provided  using  the  INSTRON  equipment  and  then  corresponding 
simulation  was  carried  out.  The  results  of  experiments  are  presented  in  Fig.  7.  Initial 
simulated  structure  is  presented  in  Fig.  8,a  (size  of  the  sample  is  1.7x0. 9  cm).  Velocity  of 
loading  plate  is  10  cm/sec.  Evolution  of  the  simulated  sample  structure  is  shown  in 
Fig.  8, a  -  c.  As  is  evident  from  presented  data  the  simulation  resiults  are  in  good 
agreement  with  those  obtained  experimentally. 


iai' 


Fig.  7.  Evolution  of  the  sample  (experimental  results). 


,  . 

1  . 

1  - 

t  .- '  i 

Fig.  8.  Evolution  of  the  sample  structure  (results  of  simulation). 
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-  Ceramics 


Based  on  Zr02  ceramic  sample  of  0.6x0,4  cm  in  size  was  simulated.  The  load  was  applied 
by  setting  up  the  constant  velocity  to  all  the  elements  of  the  upper,  layer  (piston).  It  was 
varied  m  the  range  of  50  ^  500  cm/sec.  Porosity  of  the  simulated  material  was  effectively 
taken  into  account  by  introducing  the  ruptures  of  the  links  between  some  automata,  that 
conesponds  to  micro-pores.  Macro-pores  were  simulated  by  removing  separate  automata 
off  the  sample.  The  porosity  of  the  simulated  sample  was  varied  in  a  wide  range  from  0  to 
50  %.  Analysis  of  the  obtained  results  showed  good  enough  agreement  (up  to  30%)  with 
experimental  data  [2,3],  Due  to  advantages  of  the  method  it  was  possible  to  discover  that 
the  regions  of  short  term  local  flow  may  be  formed  even  in  brittle  Zr02  (Fig.  9,a).  The  size 
of  this  region 'is  of  about  0.08  cm.  The  life-time  of  this  phenomena  is  of  about  1.0*10" 
^  sec.  After  this  time  the  region  is  smeared.  Analysis  of  the  links  between  the  automata 
showed  that  regions  of  local  flow  correspond  to  formation  of  micro-damages  (Fig.  9,b,c). 
The  calculation  showed  that  in  many  cases  the  main  crack  propagates  along  the  path 
indicated  by  the  local  flow  regions  (Fig.  9,.d). 


Fig.  9  Part  of  the  sample  in  the  region  of  local  flow;  automata  velocity  field  (a).  The  grid 
represents  the  links  of  the  automata;  b)  before  micro-damage;  c)  after  micro-damage;  d) 
after  fracture. 


80 


-  Concrete 


The  sirength  as  well  as  durabilily  of  concrete  can  be  evalualcd  by  compression  lest  of 
concrete  cyltnders^The  s.zes  of  the  standard  cylinder  are  6  inches  ,n  dLneter  atrf  12 

ref’lT  '"f  ivT/?a  experimental  slress-slram  data  arc  shown  in  Fig  10  The 

results  of  MCA  modeling  of  tlie  compression  tests  of  a  cement  paste  aggregate  and 

HemT  A  ‘‘n  presented  in  Fig.  10  have  to  be  expfained  in 

etail  As  tlie  first  step  the  compression  testing  of  cement  paste  was  simulated.  The  inter- 
automata  mteracbons  and  properties  of  a  single  cement  paste  automaton  was  fitted  to  the 
experimental  stress-strain  data  (Fig  10)  The  same  procedure  for  the  aggregate  was  used  at 
I  te  second  step.  Note,  that  the  size  of  the  both  types  of  the  automata  was  equal  to  2  0»10'^ 
cm.  At  the  last  step  the  concrete  sample  was  generated  directly  by  the  cement  paste  and 
the  aggregate  automata.  Tlie  concrete  had  an  aggregate  content  of  about  52%  The 

McHn?!  of  the  concrete  and  those  obtained  by 

MCA  modeling  (Fig.  10)  shows  a  good  agreement. 


Fig.  10.  Stress-strain  relation  for  aggregate,  cement  paste  and  concrete  solid  line 
expenmental  results;  doited  line  -  MCA  simulation  results. 


Fig.  1 1  shows  the  calculated  structure  of  the  inter-automata  linkages  and  one  of 
the  regular  kind  of  concrete  cylinder  after  compression  test  [4].  The  similar  type  of  the 
failed  pan  IS  evident.  ^ 
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2.4.  Structure  Behavior 

To  demonstrate  an  ability  of  the  developed  approach  to  describe  concrete  and 
rcinforccd-concrcte  constmctions  computer  simulation  of  beha\  ior  and  fracture  of  the  two- 
storied  slnicture  under  shcar-t>pe  deformation  imitating  crust  movement  during 
earthquake  has  been  made.  For  this  purpose  two  dimension  frame  structure  was  e.xamincd. 
Fig.  12  shows  breadboard  and  measurement  of  this  sample. 

In  the  Fig.  12  one  can  see  that  the  stationary  stopper  limits  the  movement  of  the 
upper  floor  of  the  structure.  Bottom  of  the  structure  imitates  the  foundation  lying  on  the 
earth  and  shifting  with  it  at  an  earthquake. 


Fig.  1 2  a)  breadboard  structure;  b)  concrete  structure. 

A  separate  automaton  of  the  structure  was  considered  as  the  concrete  one  So  the 
mechanical  properties  of  an  automaton  tied  with  the  same  of  the  test  concrete  cylinder.  In 
this  calculations  the  size  of  an  automaton  was  equal  to  1  cm.  Specific  characteristics  of 
cellular  automata  and  parameters  of  an  inter-element  interaction  were  measured  up  to 
with  the  parameters  of  the  high-strength  concrete  [4]. 

Displacement  of  the  crust  was  imitated  by  the  leftward  even  movement  of 
automata  of  the  construction  foundation.  The  velocity  of  the  movement  was  equal  to  1 
m/sec.  Damping  the  upper  floor  of  the  stnicture,  the  stationary  stopper  (Fig.  12)  controlled 
the  shear  type  of  the  structure  deformation. 

Fig.  13,a  shows  the  prime  stages  of  the  fracture  of  the  considered  structure  under 
the  above  mentioned  loading  conditions.  It  can  be  seen,  that  at  the  first  stage  the  cracks 
are  initiated  and  grow'  near  the  butt-joints  of  the  middle  floor  with  the  w'alls  of  the 
structure.  This  is  the  result  of  the  process,  when  the  walls  are  bonded  by  the  moving 
foundation  squeeze  "girder"  and  lend  the  wavy  contortion  to  it. 
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Further  damage  of  the  structure  proceeds  near  the  butt-joints  of  the  walls  with  the 
foundation  of  the  structure.  So,  the  structure  tends  to  tilt  as  a  whole.  At  the  same  lime 
another  crack  is  initiated  and  grows  in  the  right  wall  near  the  butt-joint  of  the  wall  with 
the  "girder".  This  is  a  result  of  the  resistance  to  bending  of  the  upper  part  of  the  right  wall 
by  the  middle  floor. 

The  developed  MCA  Method  makes  it  possible  to  analyze  the  influence 
constructive  modifications  of  the  frame  on  response  and  fracture  of  the  structure.  To 
illustrate  such  abilities  of  the  method  in  the  present  work  various  modifications  of  the 
initial  structure  (Fig.  12,b)  were  introduced. 

Fig.  12,b  shows  the  prime  stages  of  fracture  of  the  structure  with  the  thickening 
middle  floor  (here  and  further  loading  scheme  analogous  to  the  showed  one  in  the 
Fig.  12, a  was  used).  One  can  see  the  fundamental  change  of  the  fracture  pattern.  The 
upper  storey  has  been  completely  broken  down,  whereas  the  lower  floor  has  been. 
Consequently  the  broken  structure  may  be  repaired,  while  the  damaged  base  structure 
(Fig.  13, a)  cannot  be  repaired. 


b) 

Fig.  13.  Structure  of  inter-element  bonds  at  different  stages  of  fracture;  a)  control 
structure;  b)  heavy  bay  structure. 
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Initial  stmcture 
a) 


Initial  structure 
b) 


Fig.  14  Differeiil  tj-pes  of  reinforcing  of  stniclure  by  sleel  trusses:  a)  horizontal  linsses 
b)  bottom  storey  trusses;  c)  double  storey  trusses. 
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MCA  Method  allows  also  to  research  the  influence  on  behavior  and  fracture  of 
niodified  structures  by  various  modes  of  strengthening  of  the  base  structure  (Fig,  12.b). 
For  example  trussing  of  the  walls  at  the  botli  storeys  of  the  stnicture  by  steel  tmsses 
(Fig.  14, a)  causes  the  fracture  of  the  walls  and  the  middle  floor  being  between  trusses.  If 
only  the  lower  storey  of  the  structure  is  tnissed  b)’  steel  trusses  (Fig.  14.b),  then  only  the 
upper  storey  is  broken  down.  The  fracture  pattern  in  this  case  is  of  the  same  type  as  in  the 
Fig.  13,b.  Tnissing  the  both  storeys  of  the  structure  by  steel  trusses  (Fig.  14,c)  leads  to  the 
fracture  localization  only  near  the  foundation  of  the  structure.  The  walls  and  middle  floor 
of  the  structure  therewith  are  uncrippled. 

The  present  work  also  demonstrates  the  ability  of  the  MCA  Method  to  analyze 
the  influence  on  the  features  of  the  fracture  of  the  structure  by  the  type  of  material  of 
strengthening  members. 


<  M 
u  /  4 

S  (  s 


I ' 

JiS.  /  /.'xO. 

5’  )  J 


€ 


m 

ss 

XX 


i 

\  x-vj. 


,f  T . ^ 


M  <: 

t  (  $ 

--■A  X  ^ 


Initial  structure  a) 


b) 


>■< 


v>— 


X-*;  ' 

.  W’ 
XX'  ^  MX  < 

H  '' 

i  <  i 


& 


V  v.\-' 

r:;-  \  ^ 

/v,.  ?  M 

M  )  M 

'•X‘ . /• 


# 


^ 


i  i 

#,«S5®§ 

'•vV' 

'  '  V 


Fig.  15.  Different  types  of  material  of  vertical  shores:  a)  steel  shores;  b)  shores  like 
reinforced  concrete. 


A  structure  with  steel  (Fig.  15,a)  and  like  reinforced-concrete  (Fig.  15,b)  vertical 
strengthening  shores  was  modeled.  It  can  be  seen  that  the  structure  with  the  steel  iRisses 
first  damaged  in  the  lower  part  of  the  walls  (near  the  foundation  of  the  frame).  Further  the 
both  floors  of  the  structure  "burst"  in  the  butt-joints  with  the  shores. 
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Anollier  fraclure  pattern  is  seen  in  the  Fig.  15,b.  Here  first  llie  middle  floor  severs  oneself 
from  the  walls  of  the  structure,  because  the  "yielding"  misses  are  not  able  to  resist  to  waw 
contortion  of  the  "girder”.  Thereafler  the  walls  start  to  crack.  At  this  takes  place,  in 

contrast  to  the  slnictifre  in  the  Fig.  KS.a.  the  upper  floor  of  the  siructurc  is  remained 
intacl. 

The  outlined  examples  demonstrate,  tliat  the  MCA  Method  allows  to  simulate 
the  behavior  and  fracture  of  all  manner  of  the  concrete  and  reinforced-concrete 
constructions  under  various  types  of  external  loading.  The  present  approach  provides  the 
way  to  research  the  influence  of  constructive  modifications  on  vitality  of  the  construction 
as  a  whole  and  its  separate  members.  It  also  allows  to  research  the  influence  of  the  same 
modifications  on  the  features  of  fracture  process. 

The  advantages  of  MCA  Method  permit  the  use  of  this  approach  to  optimize  the 
structure  of  the  analyzed  construction  according  to  specific  requirements  arid  the  computer 
aided  expertise  of  concrete  constructions  for  natural  catastrophes. 

3.  Final  remark 


At  the  conclusion  it  should  be  pointed  out  that  because  Movable  Cellular 
Automata  Method  is  based  on  the  particle  concept  and  includes  all  the  advantages  of 
classical  cellular  automata  it  is  very  easy  to  develop  special  aided  software  for  specific 
problems  of  various  types. 
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Special  report 

PHYSICAL  MESOMECHANICS  OF  MATERIALS  AND  ITS  IMPACT  ON 

SHOCK  CHEMISTRY 

Panin  V.E,,  Makarov  P.V.,  Smolin  LY. 

Institute  of  Strength  Physics  and  Materials  Science,  SB  HAS,  Tomsk,  Russia 

Abstract.  The  fundaiiienlal  difference  of  the  physical  mcsomechanic  approach 
from  traditional  description  of  deformation  and  fracture  of  materials  consists  in  an  explicit 
form  consideration  of  internal  structure  heterogeneities  of  materials  on  the  mesoscale 
level.  Experiments  and  theory  demonstrate  that  just  these  heterogeneities  -  grains 
boundaries,  inclusions  and  etc.  -  cause  essentially  heterogeneous  fields  of  stresses  and 
deformations  in  mesovolumes  of  material  under  loading.  Just  mesolevel  of  consideration 
allows  to  clear  up  mechanisms  of  plastic  deformation  and  fracture  and  to  understand 
princij)les  of  physical  and  chemical  transformations  related  to  deformation.  As  appendix 
of  the  developed  approach  to  shock  chemistry  it  is  considered;  1)  processes  of  localization 
of  defonnation  on  the  mesolevel  and  formation  of  the  centers  of  reaction;  2)  development 
of  cumulative  processes  and  formation  of  stable  vortexes  driven  in  shock  wave  front. 

1.  Basic  principles 

The  concept  and  method  of  our  approach  are  based  on  the  ideas  of  physical 
mesomechanics  [1,2].  The  fundamental  difference  of  the  physical  meso-mechanics 
approach  from  traditional  description  of  deformation  and  fracture  process  in  materials 
consists  in  the  following; 

a)  A  material  is  considered  as  hierarchically  organized  system  of  structural  levels  of 
different  scales;  micro,  meso  and  macro. 

b)  Heterogeneities  of  internal  structure  of  materials  on  the  meso-scale  level  are  considered 
in  the  explicit  form. 

On  macroscopic  level  -  the  level  of  the  averaged  description  in  .the  model 
contributions  from  micro-  and  mesolevels  to  shear  strength  are  taken  into  account.  So  a 
macropartical  is  the  representative  mesovolume. 

On  microscopic  level  heterogeneities  of  defect  structure  are  considered  by  using 
of  defect  distribution  function  and  dislocation  kinetics  for  plastic  strain  rate  is  suggested. 
So  deformation  of  a  material  on  microlevel  is  described  in  terms  of  continuum  theor}^  of 
defects.  Taking  into  account  processes  of  nucleation  and  multiplication  of  defects  in  the 
elastic  precursors  and  fronts  of  strong  and  weak  shock  waves  kinetic  constitutive 
equations  arc  constmcled  for  dislocation  continuum. 

The  approach  allows  to  follow  the  microstmeture  evolution.  Dependencies  of 
factors  of  viscosity  of  shock  loaded  metals  are  investigated.  Evolution  of  shear  strength  in 
metals  (its  temporary  decreasing  in  shock  wave  front  and  recover)'  behind  it)  is  described 
as  process  of  localization  of  deformation  in  shear  planes  of  microscopic  scale 
(microbands).  It  is  caused  by  their  local  heating  up  and  subsequent  cooling  down. 


On  the  mesoscopic  level  heterogeneities  of  internal  mesostructure  are  considered 
in  explicit  form.  Experiments  [3]  and  theory  [4,5]  demonstrate  that  grains  boundaries, 
inclusions  and  etc.  cause  essentially  heterogeneous  stress  and  deformation  fields  in  meso- 
volumes  of  material  under  loading.  It  results  in  development  of  shear  bands  formation  and 
other  ty  pes  of  deformation  localization.  Just  mesolevel  of  consideration  allows  to  clear  up 
mechanisms  of  plastic  deformation  and  fracture  and  to  understand  principles  of  physical 
and  chemical  transfonnations  related  to  deformation.  Just  in  volumes  of  mesoscale  level 
in  materials  under  loading  areas  of  localized  deformation  are  formed  depending  on  initial 
heterogeneities  and  conditions  of  loading,  self-organizing  processes  develop.  Rotary 
plasticity  is  developed  and  material  is  found  hardly  fragmented  and  hashed. 

But  on  the  macroscale  level  of  averaged  description  the  information  connected 
with  fluctuations  and  localized  defonuation  on  the  mesolevel  is  completely  lost.  Though 
just  fluctuations  determine  further  development  of  related  to  each  other  processes  of 
chemical  and  physical  transformations  and  deformation  and  fracture. 

2.  Results  of  calculations 

A  map  of  a  polycrystalline  A1  mesovolume  is  shown  in  fig.  1.  The  results  of 
numerical  modeling  of  shock  wave  propagation  in  this  mesovolume  are  shown  in  fig.  2-4. 
On  mesolevel  the  front  of  a  shock  wave  fonns  mesoscopic  shear  bands  which  cover  from 
unit  grains  up  to  several  grains  (fig.  3).  The  process  of  formation  of  fragments  is  finished 
in  waves  of  unloading.  Opportunity  of  development  of  microcumulative  processes  and 
formation  of  microvortexes  appears.  Fg.  6  shows  derivation  of  velocities  from  average 
values  in  the  rows  of  computational  grid  perpendicular  to  the  direction  of  shock  wave 
propagation.  Stresses  corresponding  to  this  derivations  of  velocities 
(j  =  pclsu  -  0.1  ^  0.3  GPa  •  These  microvortexes  can  move  in  front  of  a  shock  wave  under 

certain  conditions.  It  should  be  mentioned  that  strain  localization  patterns  under  quasi¬ 
static  loading  in  the  same  mesovolume  is  qualitively  different  and  one  of  them  is  shown  in 
fig.  7. 


Conclusions 

Localization  of  deformation  is  caused  by  structural  elements  of  meso  scale  levels. 
They  are  tlie  grains  in  the  offered  example.  Therefore  the  length  of  the  meso  shear  bands 
is  in  order  of  grain  size.  Shear  bands  are  formed  in  the  direction  of  main  shear  stresses.  In 
some  cases  joints  of  these  shear  bands  form  bands  of  localized  defonuation  of  greater 
scale. 

The  rotational  mode  of  deformation  is  developed.  In  the  presented  example 
rotation  of  individual  grains  form  the  elementary  rotational  cells.  Rotations  of  grains 
cause  shear  bands  generation  near-by  grain  boundaries  and  grain  joints. 

So  a  new  approach  for  description  of  shock  induced  chemical  processing 
determined  by  heterogeneities  of  internal  structure  on  the  micro  and  meso  levels  is 
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presented.  It  includes  numerical  mmodeling  of  the  meso  shear  bands  formation  in  metals 
under  shock  wave  loading  and  its  impact  on  shock  induced  chemical  processing. 
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X 

Fig.  1.  A  map  of  mesovulume  of 
polycrystalline  Al. 


Fig.  2  Distribution  of  stresses  Oyy  in  the 

mesovolume  along  the  lines  parallel  to  the 
direction  of  shock  wave  propagation. 


Fig.  3.  Contour  map  of  plastic 
deformation  and  meso  shear  bands 
formation  in  the  mesovolume  of 
material  under  shock  wave 
loading. 


Fig.  4.  Distribution  of  plastic  defor¬ 
mation  in  the  mesovolume  along  the 
lines  parallel  to  the  direction  of  shock 
wave  propagation. 
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Special  reporl: 

INVESTIGATIONS  AT  MULTI-USER  ORIENTED  TsNIIMASH 
LARGE-SCALE  BALLISTIC  FACILITY 
P.V.Kryukov 

(CSRI  of  Mechanical  Engineering  Problem  .  Moscow,  Russia) 

*  OBJECTIVES  OF  THE  PROJECT  AND  THE  CURRENT  STATE  OF  THE 
ART 

The  objective  of  the  project  is  to  provide  to  various  interested  organizations  and 
individual  scientists  a  free  access  to  the  TsNIIMASH  large-scale  ballistic  facility  and 
support  them  in  conducting  their  own  research  at  it.  The  facility  has  unique 
characteristics  and  for  many  years  was  not  available  for  both  foreign  and  domestic 
civilian  investigators  being  solely  used  for  development  of  strategic  weapons.  It  is 
planned  to  use  donations  of  international  scientific  foundations  operating  at  the 
Russian  territory  and  money  allocated  from  the  Russian  National  Budget  to  prepare 
and  conduct  ballistic  experiments  within  the  project,  i.e.  to  provide  the  facility  to  the 
international  community'  of  investigators  free  of  charge. 

The  TsNIIMASH  large-scale  ballistic  facility  is  the  largest  and  the  most  complicated 
test-  bench  in  the  world  amongst  its  kind.  Each  experiment  has  a  capacity  to 
simultaneously  load  various  samples  with  different  content  which  opens  an 
opportunity-  to  exploit  the  facility  in  multi-user  oriented  mode  and  invite  a  great 
number  of  research  organizations  and  prominent  individuals  from  many  countries. 
Unique  sets  of  launch  velocities  and  loading  parameters  (pressure,  loading  cycle 
duration,  shear,  and  temperature)  can  be  achieved  at  the  facility;  no  other  facility  or 
method  can  provide  comparable  sets  of  loading  values. 

A  solid  body,  when  compressed  to  high  pressures  at  the  presence  of  certain  activating 
factors,  such  as  temperature  and  shear,  can  experience  a  vast  spectrum  of  phase  and 
other  transformations.  This  provides  an  investigator  opportunities  to  research: 

*  Processes  in  numerous  solid  phase  reactions  that  form  meta-stable  polycrystalline 
states  with  either  amorphous  or  nano-phase  structures  which  can  possess  diverse 
valuable  properties  (superconductivity  at  high  critical  parameters;  unique  magnetic, 
mechanical  or  catalytic  characteristics); 

*  Formation  of  compacts  of  powders  of  various  materials,  including  diamond 
powders  at  densities  close  to  the  x-ray  density; 

*  Processes  involving  different  polymorphous  phase  transformations,  e.g.  in  carbon, 
boron  nitride,  etc. 

In  the  framework  of  the  project  investigations  it  is  foreseen  that  advanced  materials 
with  unique  physical  and  chemical  properties  will  be  synthesized  and  subsequently 
applied  in  various  industries  on  a  large  scale. 

The  following  additional  directions  of  aeroballistic  experiments  were  selected  for  the 
TsNIIMASH  facility': 

*  Investigation  of  destruction  dynamics  of  massive  body  colliding  with  a  compact 
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inipaclor: 

*  Investigation  of  combustion  process  in  a  fuel  injected  into  a  hypersonic  air  flow; 

*  Investigation  of  scram  acceleration  of  a  Inpersonic  projectile  moving  in  the 
bore  filled  with  an  oxygen/hydrogen  mix: 

*  Investigation  of  various  methods  of  gasdynainic  acceleration  of  projectiles 
to  velocities  over  10  km/s. 

•  INVESTIGATIONS  ON  SYNTHESIS  OF  NEW  MATERIALS 

It  is  proposed  within  the  project  to  conduct  dynamic  experiments  at  the 
TsNIIMASH  large-scale  ballistic  facility  and  perform  loading  ampules 
containing  various  samples.  These  ampules  with  the  samples  are  positioned  upon 
the  target  which  is  then  installed  at  the  end  of  the  ballistic  bore  of  the  facility.  The 
dynamic  loading  of  the  samples  is  planned  to  attain  as  a  result  of  a  high-speed 
impact  of  a  30-kg  steel  disc  (impactor),  having  diameter  0.5  m  atid  20  mm  thick,  at 
a  velocity  2.5  km/s. 

In  case  that  the  ballistic  bore  is  air  evacuated,  the’  impactor  strikes  the  target 
yielding  the  maximum  pressure  around  60  Gpa  during  8  ms.  If  the  bore  is 
preliminary  filled  with  gaseous  hydrogen  tlien  the  shock  is  dampened; 
hydrogen  also  reduces  the  pressure  slope,  makes  it  smoother  as  well  as  helps 
homogenous  pressure  distribution  upon  the  target  area.  Variation  of  initial 
pressure  of  tlie  hydrogen  inside  the  ballistic  bore  and  tlie  impactor's  velocity 
provides  control  of  loading  conditions  in  the  vast  range  of  loading  pulse  length  and 
maximum  pressure  values.  At  longer  loading  pulse  and  the  pressure  higher  than 
the  yield  limit  of  the  target  material,  the  target  is  being  plastically  deformed  under 
radial  unloading  waves.  Hence,  to  preserve  the  central  partof  the  target  0.3  m  in 
diameter  from  destruction,  the  loading  pulse  length  should  be  not  greater  than 
40  ms.  Due  to  the  sound  unloading  waves  in  tlje  dampened  impact  the  profile  of 
the  pressure  in  the  bulk  of  the  samples  is  fairly  homogenous  thus  practically 
excluding  their  shock-wave  healing.  A  long  30  ms  pressure  slope,  which  is 
formed  when  the  impactor  moves  along  the  ballistic  bore,  helps  get  rid  of  cavities  in 
powder  materials  at  a  moderate  speed  without  their  collapse.  As  the  hydrogen 
adiabatically  expands  and  pushes  the  impactor  in  the  reverse  direction,  back 
from  the  target,  the  samples  in  ampules  are  being  unloaded.  The  layer  of  the 
dampening  gas  prevents  formation  of  tensile  stresses  during  samples'  unloading 
wliich  othenvise  could  have  led  to  mechanical  damage  of  the  samples.  The  second 
loading  cycle,  which  otherw'ise  could  have  been  caused  by  the  second  impact  of  the 
disc  impactor,  as  well  as  contact  of  the  samples  with  the  high  temperature  gas  jets 
are  prevented  because  a  compressed  air  is  injected  in  the  gap  between  the  recoiling 
impactor  and  the  target. 
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Parameters  Attained  on  Expereinental  Balistic  Facility 


Ballistic  bore  caliber 

0.5  m 

Ballistic  bore  length 

60  m 

Full  length  of  the  facility 

200  m 

Detonating  mix  formula 

2H2+O: 

Pressure  beyond  the  Detonation  Front 

up  to  300  atm 

Energy  yield  at  mix  combustion 

up  to  1.2  GJ 

Mass  of  Discus  Projectile/Striker 

9  kg 

Max.Terminal  Velocity  of  Discus  Projectile 

up  to  3500  m/s 

Max.Kinetic  Energy  of  Discus  Projectile 

up  to  55  MJ 

Discus  Angular  Rate 

250  rps 

Discus  Rotation  Axis/Bore  Axis  Misaligment 

better  than  30” 

Maximum  Magnetizition  of  Discus  Material 

1.4  Tesla 

High  Speed  Throwing:  Experimental  Results 

(V=2  km/s) 


Discus  projectile/striker 


Cast  iron  target 


Special  report 

CALCULATION  OF  THE  SHOCK  HUGONIOT  FOR  METALS  AND  ALLOYS. 


VM. Kuznetsov,  G.E.Rudenskii,  R.LKadyrov,  P.P.Kaminshi 
Institute  of  Strength  Physics  and  Material  Science,  SB  RAS. 

1.  Introduction 

The  problem  of  shock  Hugoniot  calculation  still  remains  to  be  a  topical  one  in 
spite  of  the  considerable  progress  achieved  in  the  field  of  theoretical  investigations  into 
solids'  properties  under  conditions  of  shock-wave  loading.  The  reason  behind  this  lies  first 
\A  ith  a  procedure  for  calculation  of  adiabatic  cun^es  where  Debye-Mie  Gruneisen  approach 
is  used,  and  second,  there  is  an  absence  of  the  reliable  enough  models  for  binding  forces 
in  solids  that  may  be  applied  not  only  to  pure  metals  but  to  their  alloys  as  well.  There  are 
various  methods  for  calculating  the  thermodynamic  properties  in  case  of  pure  metals  such 
as  pseudopotential  method  [1],  effecth  e  medium  method  [2],  and  embedded-atom  method 
[3],  nevertheless,  such  a  problem  has  not  yet  been  solved  in  case  of  transient  metals.  This 
may  be  due  to  the  specificity  of  their  atoms'  electron  structure  and  taking  it  into  account 
correctly  requires  more  complex  models  of  the  atomic  interaction  to  be  employed. 

In  the  present  paper  the  binding  forces  model  is  proposed  for  metals  and  alloys. 
This  model  is  based  on  the  formalism  of  electron  density  functional  and  allows  to  take 
into  account  many-body  interaction  effects.  The  model  parameters  are  determined  from 
the  consistency  between  the  theoretical  and  experimental  properties  of  pure  metals.  No 
extra  parameters  are  introduced  in  relation  to  alloys.  A  distinctive  feature  of  the  developed 
method  is  its  relative  simplicity’  of  calculation,  the  absence  of  externally  dependent 
parameters  what  alIo\^s  one  to  emplo\  thus  obtained  potentials  in  studying  the  kinetics  of 
processes  applying  the  molecular  dynamics  simulation. 

The  de^  eloped  method  has  been  applied  in  this  work  for  calculation  of  shock 
Hugoniot  in  Ni,  A1  metals  as  well  as  in  NiAl  (B2-structure)  and  Ni3Al  (LI 2-structure) 
alloys.  The  calculation  of  the  shock  Hugoniot  both  on  pure  metals  and  alloys  has  been 
conducted  ith  the  use  of  the  quasiharmonic  approximation  without  resort  to  the  Debye- 
Mie  Gruneisen  model. 

A  calculation  procedure  has  been  proposed  to  obtain  the  result  that  can  be 
applied  to  investigate  the  phase  transitions  in  alloys  and  to  develop  a  plan  of  shock-wave 
experiment. 

2.  Electron  density  functional  method 

MEDF  procedure  is  based  on  the  theory  of  electron  density-  functional  where  the 
full  energ\'  of  a  solid  body  is  represented  by  functional  E[p]  of  the  electron  density-  p(r). 
The  exact  form  of  E[pl  expression  is  not  known  in  case  of  solids. 

To  construct  a  model  electron  density  functional,  the  configurational  expansion 
has  been  applied  to  the  many-electron  function  in  terms  of  one-electron  determinants 
w  hich  basis  functions  were  selected  using  the  Wannier  functions'  representation.  This 
giv  es  us  an  opportunity  to  introduce  (in  terms  of  one-electron  approximation)  non-integer 
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filling  numbers  for  tlie  electron  shells  of  atoms  as  well  as  to  take  into  account  sp- 
hybridization  effects  in  transition  metals.  Such  an  approach  enables  us  to  divide  all  the 
electrons  in  the  system  into  core  ones  (internal  atom  shells  and  valence  d-electrons  in  case 
of  tlie  transition  metals)  and  nearly  free  ones  what  allows  us  to  write  well-founded  enough 
expressions  for  subsystems'  energy'.  A  criterion  for  dividing  the  electrons  into  two 
subsystems  is  a  degree  of  localization  of  the  considered  electron  shells  in  relation  to  the 
atom  nucleus.  In  this  work,  for  instance,  the  number  of  free  electrons  in  A1  has  been  taken 
equal  to  3  (configuration  3s^3p*)  and  1.4  for  Ni  (configuration  3d^  ^4s’ '')  in  accordance 
^vith  the  results  of  quantum  mechanics  calculation. 

On  such  subdivision  of  electrons  into  two  subsystems,  full  energ\'  of  a  solid  body 
mav  be  written  in  the  following  form; 

E  =  EJpJ+EJpJ+E,[p^,pJ .  (1) 

where  first  component  represents  energv’  of  free  electrons,  second  component  is  energv'  of 
core  electrons  and  third  component  relates  to  the  interaction  beUveen  the  above  mentioned 
subsystems. 

Nearly  free  electrons. 

As  it  has  been  done  according  to  the  pseudopotential  norm-conserving  method 
[1]  let  us  move  from  the  real  density  of  nearly  free  electrons  that  suffers  intense 
oscillations  in  the  vicinity'  of  atom  nucleus  due  to  it  being  orthogonal  to  core  electrons  to 
pseudo-wave  functions  that  only  fluently  decrease  there.  Due  to  passing  to  the  pseudo¬ 
electrons  that  possess  corresponding  distribution  of  their  density'  pjf)  we  now  can  write 

their  Ml  energv  using  the  nearly  uniform  electron  gas  approximation  [4],  namely: 

A  [PJ  =  [P.  ]  +  TJp.]  +  A?:,  [pj  +  C  [pj  +  ET'  [ A  7  (2) 

Here  £  and  are  successively  Coulomb,  kinetic  and  exchange  components  of 
the  valence  subsystem  energ> AT;,  is  the  Kirzhnitz’s  correction  [5]  for  nonuniformity'  of 

the  electron  density,  and  El''"  is  the  correlation  energv  in  the  form  obtained  by  Gombas 
[6]  as  a  result  of  extrapolation  of  electron  gas  correlation  energ>'  in  case  of  high  and  low 
density'. 

To  approximate  the  density'  of  valence  pseudo-electrons  in  MEDF  one  may  use 
the  approximation  of  superimposed  spherically  symmetric  atom  functions  of  density'  ; 

A  (r)  =  I C; (r-R,)^  I  A“  N }dV  =  ,  (3) 

j.a 

where  is  a  number  of  free  electrons  in  atom  of  the  type  a. 

The  second  order  smallness  expansion  of  nearly  free  electrons'  energv 
components  (2)  as  a  function  the  free  electrons'  mean  density'  allows  us  to  present  E,.  as 
follows. 


98 


(4) 


EJpJ^EJpJ  +  ^YYCrcf  v:^  (p^ 

ccp  i,j 

e  where  energy  Eo  as  well  as  appearing  in  (4)  effective  pair  potential  are  dependent  on 

the  mean  density  of  valence  electrons  (on  crystal's  volume).  The  value  of  equals 

to  1  or  0  depending  upon  whether  an  i-lattice  point  is  taken  by  the  atom  of  the  type  a 
or  not. 

Core  electrons. 

The  energ>'  of  core-core  interaction  is  given  using  the  Kim-Gordon 
approximation  [7]  as  a  sum  of  effective  pair  short-range  potentials  Vap(R,j),  which  include 
(like  it  is  done  with  the  tight-binding  method)  Coulomb,  kinetic,  exchange  and  correlation 
components: 

EJpJ^^YYCrcX/R^j).  (5) 

ap  ij 

The  energy  of  isolated  cores  is  taken  as  an  initial  point  in  fiill  energy  calculation.  The  Ec 
contribution  to  the  binding  energ>’^  is  zero  in  case  that  densities  of  core  electron  of 
neighboring  atoms  are  not  overlapped. 

Pseudopotential 

The  energ>'  of  core-valence  interaction  in  case  of  alloy  take  the  form: 

E^  =  l.q\drW,(\r-Rj\,pPlpJr).  (6) 

j,a 

Here  Wa  is  a  local  norm-conserving  interaction  pseudopotential  for  atoms  of  the  t>pe  a  .  It 
has  an  effect  on  free  electron  density  pv  moving  to  r.  Let  us  note  that  the  pseudopotentiai 
is  explicitly  dependent  on  free  electron  density  pv ,  whereas  the  density  dependence  of  the 
core  electrons  is  a  pure  parametrical  one  because  of  we  have  used  frozen  cores 
approximation. 

The  full  variation  problem  in  case  of  s>^stem  of  particles  is  in  finding  such  type  of 
pv  using  which  minimizes  energy  E  under  given  external  conditions.  This  is  equivalent  to 
the  solution  of  the  problem  5E=0.  To  avoid  difficulties  stemming  from  the  necessity  to 
find  the  form  p,.  for  each  new  atomic  configuration,  the  pseudopotential  is  constructed 
in  such  a  manner  that  a  fiill  energy  minimum  condition  for  any  type  of  configuration 

I  Rj  I  would  be  true  for  one  and  the  same  fixed  type  of  nearly  free  electron  density 
function  (3). 

Being  the  result  of  expansion  of  the  core-valence  interaction  up  to  second  order 
of  smallness  and  in  the  vicinity  of  (equilibrium  mean  nearly  free  electrons'  density 
of  the  metal  of  a  type  )  full  energy  may  be  written  in  the  form: 
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(7) 


E  =  EJpJ  +  \Y.Y.Cr  (p, ,  R,  )  + 

aP  i,j 

uPy 

where  the  second  term  includes  in  itself  pair  energy  components  of  both  core  and  valence 
subsystems  together  with  their  interactions  and  the  last  term  is  a  three-particle  addition  to 
the  energ}'  of  the  core-valence  interaction. 

To  determine  pseudopotential's  parameters  we  resorted  to  the  following  known 
from  the  experiments  properties  of  pure  metals:  equilibrium  volume,  bulk  modulus,  shear 
modulus  C44,  vacancy  formation  energt',  pressure  at  20%  of  compression  and  tension  at 
zero  isotherm  being  approximated  with  the  use  of  the  universal  equation  of  state  [8], 

To  minimize  the  fiill  energ\',  the  nearly  free  electron  density  function  pv  is 
constructed,  according  to  (3),  as  a  superposition  of  the  atomic  density  functions  of  alloy 
components  being  approximated  by  Slater's  functions.  The  Slater's  functions'  parameters 
were  selected  in  such  a  manner  that  to  provide  the  best  coincidence  with  tabulated  by 
Herman  and  Skillman  [9]  electron  density  distribution  for  separate  atoms. 


Alloys.  • 

The  full  energ>'  of  an  alloy  is  dependent  not  only  on  its  atomic  volume  but  also  on 
its  crv'stal  structure  parameters.  Because  of  this,  to  determine  basic  state  thermodynamic 
characteristics  of  NiAl  and  NisAl  alloys  there  was  carried  out  a  procedure  for  energ\^ 
minimization  o\  eT  the  atomic  volume  and  ciy^stal  structure  parameters. 

The  basic  state  of  NisAl  is  Lty  structure  type  with  the  atomic  volume  vo=76.12 
a.u.  A  competitive  phase  is  f  c.c.  DO22  structure  type.  Moreover,  we  considered  also  its 
f  c.c.  lattice  —  base  structural  modification  DO3,  The  calculations  of  made  by  us  for  NisAl 
formation  enthalpy  of  each  considered  structural  modification  showed  that  LI 2  phase  is 
stable  at  the  equilibrium  volume  Vo=76.38  a.u. 

Basic  state  of  NiAl  alloy  is  realized  as  B2  structure  type  with  vo-82.28  a.u.  The 
calculation  of  NiAl  formation  enthalpy  in  case  of  competitive  structure  types  b.c.c.  B32 
and  fc.c.  Llo  have  supported  energy  stableness  of  B2  phase  with  atomic  volume  of 
82.26  a.u. 

3.  Elastic  moduli  and  phonon  spectra  calculations. 

To  verity'  the  applicability'  of  MEDF  for  description  of  thermodynamic  properties 
of  the  considered  metal  and  alloys,  phonon  spectra  and  Fuks's  second  order  elastic  moduli 
(the  so-called  static  moduli)  w'ere  calculated  for  Ni.  Al,  NiAl  u  NtyAl.  The  calculation 
results  for  the  bulk  modulus  B  and  shear  moduli  C44  and  C*  (at  T=0  and  P=0)  are  given 
below  in  Table  1  together  with  the  corresponding  values  obtained  from  experiments. 
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Table.  1.  Elastic  moduli  of  Ni,  Al,  NiAl  and  Ni^Al.  Experimentally  obtained  values  were 


borrowed  from  flOJ.  All  values  are  given  in  GPa. 


B 

Bexp 

C44 

C44exp 

C’ 

C’exp 

Al 

79 

79 

56 

31 

21 

26 

Ni 

187 

188 

118 

132 

30 

55 

NiAl 

156 

166 

101 

112 

23 

35 

NiaAl 

187 

175 

108 

129 

30 

40 

The  reduced  value  of  C  modulus  in  all  considered  systems  is  due  to  usiiig  approximation 
of  spherically  symmetrical  atomic  functions  of  density  for  core  and  valence  electrons. 

Since  contribution  from  the  thermal  oscillations  was  taken  into  account  in 
calculating  shock  Hugoniot  by  a  quasiharmonic  model  where  the  wave  vector  k 
dependence  of  oscillations'  frequency  is  used,  we  calculated  the  phonon  dispersion 
dependencies  o)j(k)  for  all  considered  above  metals  and  alloys.  Tlie  phonon  spectra 
calculation  results  for  NiAl  and  NijAl  are  illustrated  in  Fig.l  in  comparison  with 
experimental  data. 

In  case  of  pure  metals  the  deviation  between  theoretical  and  experimentally 
obtained  dispersion  curves  does  not  exceed  10%. 


4.  Calculation  of  equation  of  state  and  shock  Hugoniot. 

The  Hugoniot  equation  takes  the  following  form  in  case  of  that  the  substance  is  at 
rest  in  front  of  the  shock  wave: 

P„ (v,T)  =  2(E„ (v,T)-EJ /  (V,  - V) ,  (8) 

where  Ph  is  a  Hugoniot  pressure: 

■  dv  .  (9) 

Here  Vo  is  equilibrium  volume  of  an  elementar>'  cell,  Eo  is  a  substance'  energy  before  the 
shock  wave  front ,  Eh  and  F  are  correspondingly  energy  and  free  energj^  of  the  substance 
beyond  the  shock  wave  front  calculated  per  one  atom.  In  case  of  quasiharmonic 
approximation 


fi 

F( v,T)  =  EJ  v)  +  v^-jdo}-(i)- gjco)  +  vJcJldcD- g„(o))-ln 


\-e 


J 


E( V, T)  =  E/ v)+Vp-ldco-o)  gJ(a)+  v/fda-  , 

g„f  /vjis  a  frequency  distribution  density  being  normalized  to  unit: 


(10) 


Vp  is  a  number  of  phonon  modes  (3  for  metals,  6  for  NiAl  of  the  B2  structure  and  12  for 
NisAl  of  the  LI2  structure);  Vk  is  a  volume  of  a  reciprocal  lattice. 

Mie-Gmneisen  approximation  is  common  to  be  used  for  calculation  of  Ph, 
namely,  the  contribution  from  temperature-dependent  part  of  F  to  the  pressure  is 
determined  by  the  expression  Pt=(Y  )A',  whereas  different  dependencies  such  as  for 
example  Landau-Slater,  Zubarev-Vaschenko  and  Dougdale-McDonald  formulas  are  used 
to  define  the  volume  dependence  of  the  Gruneisen  parameter.  Apart  from  Mie  -Gruneisen 
approximation,  Debye's  approximation  Et(v,T)=3kBTD(0/T)  may  be  applied  to  calculate 
temperature  dependence  of  energ>^Et(v,T),  There  is  a  number  of  various  dependencies  to 

approximate  volume  dependence  of  the  Debye  temperature  e(v),  one  of  them  is  an 
approximation  based  on  representation  of  y  via  Debye  frequency  cod;  y-51na)d/5inv  from 
where  one  can  derive  the  dependence  cod(v)  under  condition  that  the  dependence  y(v)  and 
equilibrium  value  of  G)d(vo)  are  known. 

A1  these  thermodynamic  models  have  a  ro^^^  of  the  following  drawbacks  which 
become  even  more  substantial  when  models  are  applied  to  calculation  of  alloy  properties. 

1)  Using  the  Debye  model,  it  is  unavoidable  to  know  the  volume  dependence  of 
tod.  The  choice  of  the  dependence  form  being  well-ground  and  reasonable  enough  still 
remains  to  be  to  a  great  extent  arbitrary  .  Furthermore,  to  use  the  model  it  is  necessary'  to 
know  the  value  of  the  Debye  frequency  ©0  at  \-=Vo  (or  yo=y(vo)  for  some  models). 
Although  its  value  was  reliably  established  from  indirect  experiments  for  pure  metals 
including  Ni  and  Al,  there  are  no  not  only  such  data  in  case  of  alloys  but  even  no  reliable 
evaluation  criteria  exist. 

2)  Substitution  of  the  phonon  spectra  gn(©)  with  Debye  density  distribution 
function  gd(©)  looks  like  very'  rough  approximation  for  being  applied  in  case  of  alloys 
because  there  are  isolated  from  acoustic  levels  peaks  of  oscillations'  optical  modes. 

3)  Applying  Mie-Gmneisen  approximation  also  results  in  necessity  to  use 
arbitrary'  approximations  for  the  volume  dependence  of  y  and  in  some  cases  it  is  required 
to  know  ^'alue  of  yo  which  has  never  been  determined  for  alloys. 

In  connection  with  the  aforesaid  we  ha^’e  made  a  decision  to  refuse  from 
employing  Debye-Mie-Gmneisen  approximations  and  conducted  the  calculations  of  j 
equations  of  state  on  metals  and  alloys  according  to  purely  quasiharmonic  model. 

Calculation  procedure  was  as  follows: 

—  gn((i5)  function  was  calculated  (linear  tetrahedron  numerical  integration 
method  [14]  following  by  ZB)  at  different  values  of  volume  v  taken  from  within  the 
interval  [0.5vo,  vo]  ; 

—  the  equation  (8)  was  numerically  resolved  (by  modified  Newlon  method)  in 
relation  to  T  at  each  value  of  volume  v,  in  so  doing  pressure  Ph  was  calculated  by 
numerical  differentiating  of  free  energy'  over  the  volume.  Energy  E  in  itself  as  well  as  free 
energv’  F  were  calculated  according  to  expressions  (10)  (integration  by  method  of  trapeze  ). 

It  is  worth  to  note  that  this  procedure  has  some  methodological  details  relating  to 
the  necessity  for  self-consistent  choice  of  a  calculation  accuracy:  a)  accuracy  of  g,,!©) 
function  which  is  defined  simultaneously  by  the  number  of  division  tetrahedrons  ZB  and 
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step  of  co;  b)  accuracy  of  pressure  Ph  which  is  a  function  of  step  of  v;  c)  accuracy  of 
Hugoniot  equation  solution  (8)  which  is  influenced  by  the  step  of  T.  When  the  choice  of 
these  value  is  done  in  a  non-self-consistent  manner,  there  could  be  observed  such  Ph(v) 
behavior  (oscillating,  spasmodic)  that  has  no  physical  meaning. 

The  results  on  calculating  shock  Hugoniot  in  Ni,  Al,  NiAl  and  NiaAl  are  shown 
in  Fig.2  and  Fig.3,  Comparing  them  to  experimental  data  for  pure  metals  one  may  see  a 
substantial  divergence  at  higher  compression  values  (x<0.65).  This  may  be  related  to  the 
frozen  core  approximation  since  Kim-Gordon  approximation  [7]  fails  in  attempt  to 
describe  such  high  compression.  As  one  would  expect,  the  results  of  the  Hugoniot 
temperature  calculation  for  pure  Al  and  Ni  are  in  good  agreement  with  half-empirical 
evaluations  of  [15-16]  shown  in  Fig.2. 

It  is  well  known  that  the  electron  conductivity  contribution  to  energ>'  and 
pressure  become  a  prominent  one  at  high  temperatures.  Free  energy  of  the  thermal 
electrons  may  be  evaluated  with  the  use  of  the  half-empirical  formula  [14]: 

,  where  Pois  an  electron  analogue  of  Gruneisen  parameter.  For 

T-^oo  Ye  tends  to  2/3.  Fortunately,  electron  contribution  to  the  pressure  and  temperature  at 
shock  Hugoniot  is  substantial  only  at  high  compression  values  (x<0.65). 

To  enable  a  comparison  with  used  in  [17]  expressions  for  y(v),  we  have 
calculated  mean  Gruneisen  parameter  Y=:(PtV)/Et  at  the  evety  point  belonging  to  a  shock 
Hugoniot.  As  it  may  be  seen  from  Fig.2.  these  results  are  well  agree  to  each  other  within 
the  compression  value  field  x>0.65. 

We  have  not  come  across  the  experimental  data  for  principal  Hugoniot  obsen^ed 
on  NiAl  and  NiaAl  alloys.  The  results  of  our  calculations  (Fig.3.)  demonstrate  more 
complicate  dependence  Ph(v)  as  compared  to  what  is  possible  to  obtain  through  use  of 
well-know  rule  for  averaging  the  Hugoniot  pressure  of  pure  metals  [13].  Unusual  behavior 
of  the  Gruneisen  parameter  in  alloys  may  be  related  to  the  specific  character  of  their 
phonon  spectra.  As  could  be  shown,  the  Gruneisen  parameter  Ymax^^lncOmax/Slnv  calculated 
from  the  maximum  frequenc>^  ©max  of  phonon  spectrum  would  be  closer  to  yai  and 
would  have  a  similar  linear  character.  Analogously  calculated  but  with  the  use  of  minimal 
frequency  (but  not  zero)  Gruneisen  parameter  Ymm  would  demonstrate  intense  oscillations 
that  would  be  similar  to  those  shown  in  Fig.2  as  average  y  for  alloys.  The  same  could  be 
said  about  small  deflection  of  Ni  y(v)  dependence.  The  result  would  be  close  to  the  exact 
value  of  Y  if  one  would  make  calculations  taking  into  account  a  medium  frequency  of  a 
spectrum.  These  regularities  will  result  in  appearing  almost  linear  segment  at  pressure  Ph 
curv'e  for  alloys  within  the  range  of  x~0.6 — 0.65  since  the  volume  dependence  for  optical 
(maximal)  differs  from  that  of  for  acoustic  (minimal)  frequencies. 

5.  Conclusion 

The  carried  out  calculations  of  shock  adiabatic  curves  on  NiAl  and  Ni3Al 
revealed  the  unusual  behavior  of  the  Gruneisen  parameter  which  may  be  related  to  the 
specificity  of  volume  dependence  of  crystal  lattice  atom  oscillation  frequency. 
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A  direct  logical  extension  of  this  work  is  an  investigation  into  structural  phase 
transitions  under  shock  loading,  calculation  of  isotropy  curves  in  unloading  of  alloys. 

The  possibilities  of  MEDF  are  not  limited  only  by  metallic  bond  systems:  the 
interaction  potentials  obtained  within  the  frames  of  MEDF  can  be  applied  for 
calculations  of  thermodynamic  properties  in  ceramic  materials. 
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Fig.l  Calculated  phonon  dispersion  curves  of  NiAl  (B2)  and  Ni^Al 
(LI2)  (solid  lines  -  longitudinal,  dashed  lines  -  transverse  modes). 
Experimental  data  (points)  from  [11]  -  for  NiAl  and  [12]  -  for 
NisAL 


Fig.2  Calculated  shock  Hugoniot  of  Ni  and  Al  {Ph  -  principal 
Hugoniot,  -  zero  isotherm).  Experimental  data  (points)  from 
[17].  Hugoniot  temperature  T  calculations  (triangles)  from  [15]  - 
for  Ni  and  [16]  -  for  Al. 
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for  Al,  dot  and  dash  -  for  Ni. 
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